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Requirements  for  the  Degree  of  Doctor  of  Philosophy 

DRIFT  FIELDS  AND  BOUNDARY  CONDITIONS  AT  GRAIN  BOUNDARIES 
IN  BIPOLAR  P/N  JUNCTION  SILICON  DEVICES 

By 

Yueh-Hsin  Chen 
August  1984 

Chairman:  Fredrik  A.  Lindholm 

Major  Department:  Electrical  Engineering 

In  this  study  we  present  an  analytical  treatment  of  drift-field 
solar  cells  and  more  generally  of  bipolar  p/n  junction  devices 
containing  a built-in  drift  field  in  which  for  the  first  time  the 
analysis  includes  the  dependence  of  mi  nority-carrier  lifetime  on 
position  in  accordance  with  a power-law  dependence  on  the  shallow 
impurity  concentration.  Design  implications  are  discussed 

quantitatively  for  a thin  drift-field  Si  p/n(x)  junction  solar  cell. 

Experimental  data  on  the  performance  of  non-optimal  drift-field 
solar  cells  and  discussions  on  determination  of  critical  parameters  such 
as  minority-carrier  lifetime  by  several  different  methods  are  also 
included. 
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We  also  make  efforts  toward  solving  the  current-voltage 
characteristics  of  a semiconductor  filament  in  the  steady  state  excited 
by  a time-independent  and  space-independent  carrier  density  injected  at 
the  front  surface.  First,  we  assume  a model  having  infinite  surface 
recombination  velocity  (zero  excess  carrier  density)  on  the  lateral 
surfaces,  a model  used  by  previous  workers  for  both  transient  and 
steady-state  conditions.  For  the  steady  state,  we  demonstrate  that  the 
exact  solution  for  total  current  crossing  the  injecting  face  diverges. 
Use  of  partial  sums  derived  from  truncating  the  Fourier  or  Bessel - 
function  series  of  the  periodic  extension  in  space  of  the  constant 
carrier  density  on  the  injecting  face  yields  closed-form  analytic 
solutions  for  the  current-vol tage  characteristic.  Second,  we  assume  the 
surface  recombination  velocity  is  finite  (106  cm/s)  on  the  lateral 
surfaces.  The  assumption  leads  to  an  infinite-series  solution.  We 
discuss  the  possibility  of  convergence  of  the  infinite  series  and 
suggest  the  truncation  of  it  to  obtain  the  current-vol  tage 
characteristic. 
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CHAPTER  1 
INTRODUCTION 


The  purpose  of  this  study  is  to  explore  theoretically,  and  to  some 
extent  experimentally,  aspects  of  silicon  p/n  junction  devices  relevant 
to  their  use  as  diodes,  solar  cells  and  transistors.  Solar-cell 
performance  provides  the  main  motives  for  much  of  the  theoretical  work 
done,  and  the  experimental  interpretation  reported  is  for  the  measured 
performance  of  Si  solar  cells. 

We  emphasize  devices  in  which  a gradient  of  chemical  potential, 
resulting  from  a gradient  of  shallow-level  dopant  concentration, 
exists.  This  gives  rise  to  an  electric  field  in  quasi -neutral  regions 
of  the  device  that  is  present  in  equilibrium  and  is  perturbed 
negligibly,  insofar  as  minority  carriers  are  concerned,  for  low- 
injection  conditions.  Kroemer  first  suggested  the  utility  of  such  a 
built-in  drift  field  in  his  proposal  of  the  drift  transistor  [1], 
Kroemer  mainly  emphasized  the  decrease  in  transit  time  resulting  from  an 
aiding  built-in  drift  field.  The  emphasis  in  this  study  rather  is  on  a 
retarding  drift  field,  both  in  transistors  to  some  extent  but  mainly  in 
solar  cells.  Wolf  first  suggested  the  drift-field  solar  cell  [2].  Sah 
and  Lindholm  [3]  have  recently  indicated  the  advantages  of  such  a 
structure,  through  exact  finite-difference  computer  solutions  using  the 
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Sah  equivalent  circuit  model  [4].  The  advantages  do  not  require  a drift 
field  across  the  entire  quasi-neutral  base  of  the  solar  cell;  rather 
they  exist  if  the  drift  field  extends  at  least  20  microns  from  the  back 
surface.  Thus  Sah  and  Lindholm  have  generalized  the  drift-field 
structure  and  connected  it  to  the  widely  used  back-surface-field  (BSF) 
structure  of  Mandelkorn  and  Lamneck,  Jr.  [5],  Sah  and  Lindholm  call 
this  generalized  drift-field  (DF)  solar  cell  an  extended  back-surface- 
field  ( EXBSF)  solar  cell. 

In  this  study  we  present  an  analytical  treatment  of  drift-field 
solar  cells  and  more  generally  of  bipolar  p/n  junction  devices 
containing  a built-in  drift  field  in  which  a correlation  for  the  first 
time  is  included  between  the  minority-carrier  lifetime  and  the  shallow 
impurity  concentration. 

Chapter  3 reports  experimental  performance  of  non-optimal  drift- 
field  solar  cells  and  discusses  determination  of  critical  parameters 
such  as  mi nority-carrier  lifetime  by  several  different  methods. 

Chapter  4 turns  to  the  study  of  polysilicon  devices  in  which  the 
grain  boundaries  are  assumed  to  have  infinite  surface  recombination 
velocity.  The  treatment  features  an  evaluation  of  the  Shockley  filament 
theory  [6]  and  demonstrates  that  earlier  numerically  based  studies  that 
used  a similar  model  [7]  could  not  have  yielded  the  results  published 
because  the  solution  for  the  model  used  diverges.  Indeed,  the  study  of 
convergence  in  Chapter  4 gave  rise  to  a new  method  of  analysis.  This 
method  alters  the  boundary  conditions  of  a model  having  infinite  surface 
recombination  velocity  at  the  grain  boundaries  so  that  convergence  is 
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assured.  The  alteration  makes  use  of  truncated  Fourier  or  Bessel  series 
realizations  of  the  periodic  extension  of  the  boundary  conditions  at  the 
p/n  junction  space-charge  region  under  forward  voltage.  Not  only  does 
the  treatment  of  Chapter  4 apply  to  the  current-voltage  characteristics 
of  junction  semiconductor  devices;  it  applies  also  to  various  analogous 
problems  in  heat  conduction  and  potential  theory. 

In  Chapter  5 appears  an  extension  of  Chapter  4,  in  which  the 
surface  recombination  velocity  at  the  grain  boundaries  is  finite.  The 
method  treated  is  that  of  Gonzalez  [8].  Chapter  5 deals  with  the 
details  of  Gonzalez's  method,  discussing  convergence. 

In  Chapter  6,  we  state  conclusions  and  recommend  future  related 

work. 


CHAPTER  2 

CURRENT-VOLTAGE  CHARACTERISTIC  FOR  BIPOLAR  P/N  JUNCTION 
DEVICES  WITH  DRIFT  FIELDS,  INCLUDING  CORELATION  BETWEEN 
CARRIER  LIFETIMES  AND  SHALLOW-IMPURITY  CONCENTRATION 

2.1  Introduction 

Most  p/n  junction  devices  contain  a position-dependent  doping 
concentration.  For  diffused  or  ion-implanted  emitter  regions  of  diodes, 
detectors,  transistors,  or  solar  cells,  for  example,  this  position 
dependence  results  from  the  fabrication  processing  used  in  forming  the 
emitter.  In  the  drift-field  junction  transistor  proposed  and  analyzed 
first  by  Kroemer,  a doping  profile  in  the  base  region  produces  a 
gradient  of  the  chemical  potential  and  an  accompanying  electric  field 
that  aids  injected  minority  carriers  across  the  base  region,  improving 
the  device  performance  [1].  In  p+/n/n+  structures,  as  another  example, 
the  n/n+  (low-high)  junction  provides  a local ized  built-in  electric 
field  that  retards  the  entrance  of  holes  into  the  n+  layer.  This 
structure,  analyzed  first  by  Hall  [9],  who  called  it  a donor  contact, 
was  first  used  to  improve  the  characteristics  of  power  rectifiers.  If 
the  n+  layer  is  thick  compared  with  the  diffusion  length  of  drift- 
modified  diffusion  length,  and  if  this  layer  has  deep-level 
recombination  center  concentrations  and  shallow-level  donor 
concentrations  low  enough  that  the  recombination  rate  is  sufficiently 
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small,  then  the  low-high  junction  in  a solar  cell  helps  prevent  the  loss 
of  photogenerated  hole-electron  pairs  by  recombination.  Thus  the 
conversion  efficiency  of  the  back-surface-field  (BSF)  solar  cell  [10] 
often  exceeds  that  of  a back-ohmic-contact  (BOC)  solar  cell.  On  the 
other  hand,  if  the  n+  layer  is  relatively  thin  and  borders  on  an  ohmic 
contact  or  has  recombination  rates  that  are  large  either  because  of 
band-bound-band  transitions  or  of  phonon-assisted  impact  Auger 
transitions,  then  a good  ohmic  contact  results  from  the  n/n+/metal 
structure.  In  the  drift-field  (DF)  solar  cell  first  analyzed  by  Wolf 
[2],  the  variation  of  donor  concentration  with  position  in  a p/n(x) 
structure  is  such  that  a retarding  electric  field  distributed  throughout 
the  base  reduces  the  base  recombination  current,  increasing  thereby  the 
power  conversion  efficiency  EFF. 

Most  previous  analytical  treatments  of  devices  having  built-in 
electric  fields  assume  a position  independent  carrier  lifetime.  Several 
exceptions  exist.  Some  workers  [11-13]  analyzed  drift-field  solar  cells 
and  related  devices  assuming  a minority-carrier  lifetime  that  is  linear 
in  position  x within  the  quasi -neutral  base  region.  This  gives  a 
general  trend  of  lifetime  decreasing  with  increasing  base  doping 
concentration,  which  is  consistent  with  gettering  of  impurities  in  the 
relatively  highly  doped  portion  of  the  base.  Although  this  particular 
x-dependence  of  the  lifetime  gives  the  correct  trend,  it  apparently  has 
no  other  phenomenological  bases.  Lindholm,  Li  and  Sah  [14]  analyzed  the 
emitter  region  of  n(x)/p  structures  utilizing  the  physical  chemistry  of 
di vacancy-donor  reactions  reviewed  by  Fuller  [15].  In  [14]  and  in  a 
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later  related  paper  [16],  a power-law  relation  (lifetime  proportional  to 
Nqq(x)  m)  resulted,  but  these  authors  provided  only  a semi -quantitati ve 
solution  of  the  relevant  hole  continuity  equation  [14]. 

In  the  present  study,  we  give  apparently  the  first  analytic 
solution  of  the  minority-carrier  continuity  equation  for  a quasi-neutral 
base  region  (assumed  n-type)  in  which  the  lifetime  is  correlated  with 
the  position  dependence  and  the  concentration  of  shallow-level 
impurities  according  to  the  physical  chemistry  reviewed  by  Fuller 
[15].  The  assumed  concentration  dependence  of  the  lifetime  is  also 
consistent  with  the  empirical  dependencies  of  Kendall  [17]  and  with  a 
recent  treatment  by  Fossum  and  Lee  [18]  containing  a review  of  measured 
maximum  lifetime  dependencies  for  n-type  Si.  It  is  also  consistent  with 
the  work  of  Sah  and  Wang  [19]  which  attributed  the  position  dependence 
of  defect  centers  near  surfaces  to  vacancy  complexes  diffused  into  Si 
from  surface  sources.  The  analytic  solution  to  be  obtained  below  has 
design  implications  for  transistors,  diodes,  and  solar  cells. 

2.2  Analysis 

We  assume  the  following  for  the  quasi -neutral  graded  n-type  region 
under  study  (the  symbols  used  have  their  usual  meaning): 


(i) 

Tp(x)  =>  K[NDD(x)rm 

(2.1) 

(ii) 

K = 2.81  x 1012  s cm-2  , for  m = 1 [18] 

(iii) 

d(Jp/q) 

0 ■ dx  + P'V*> 

(2.2) 
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(iv)  Jp(x)  = qiipP(x)Exo  - qDpM*i  (2.3) 

X 

(v)  exp[-  / (qExo(z)/kT)dz]  = Ndd(x)/Ndd(0)  . (2.4) 

Assumption  (i),  with  m=l,  is  consistent  with  [15,18,19]  for  n-type  Si 
for  Ndd  » 7 x 1015  cm"3  and  is  optimistic  (predicts  too  large  a 

1 C _ O 

lifetime)  for  <7  x 10  cm  . Assumption  (i),  with  m=2, 
corresponds  to  band-band  phonon-assisted  impact  Auger  recombination, 
which  can  become  the  dominant  process  for  Si  for  >„  7 x 10^  cm-^. 
Assumptions  (i)  and  (ii)  are  appropriate  for  Si  that  is  effectively 
gettered  [18].  Assumption  (iii)  is  the  hole  continuity  equation  for 
low-injection  conditions,  for  which  we  will  provide  new  analytic 
solutions.  Assumption  (iv)  expresses  minority  hole  current  density  as  a 
sum  of  drift  and  diffusion  components.  All  x-dependencies  are 
explicitly  indicated.  The  mobility  ^ and  diffusivity  Dp  vary 

sub  1 i nearly  with  Nqq  whereas  Tp  varies  linearly  or  superl i nearly . We 
assume  therefore  that  appropriate  average  values,  jip  and  Dp,  to  be 
defined  later,  will  yield  solutions  of  (2.2)  that  are  not  misleadingly 
inaccurate.  Assumption  (v)  gives  the  built-in  electric  field  Exo  for 
full  ionization  of  the  shallow  donor  atoms  for  Boltzmann  statistics;  EYO 
is  nearly  space  independent  subject  to  the  assumption  of  an 
approximately  exponential  doping  profile. 

Much  evidence  suggests  m=l  in  (2.1)  [17,18].  We  shall  use  m as  a 
general  parameter  in  solving  (2.1)-( 2.4) , but  will  use  m=l  in  discussing 
device  design. 
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If  Ndd(x)  = NDD(0)exp(Gx) , (2.4)  becomes 


Exo  - -kTG/q  . 


Combining  (2.1),  (2.2),  (2.3),  and  (2.5),  we  have 


d^P(x)  _ qExo  dP(x)  LI1PP 
dx2  kT  dx  KD 


CNnn(0)3 


m 


[exp(mqExo/kT)]P(x)  = 0 


Let  P(x)  = h(y)  and  y = exp(mqEx0/2kT) 
as  suggested  in  [20].  Then  (2.6)  becomes 


2 2C2 
m q E 


4k2T2 


2 9 m2  2r2  2C2  rM  ,nnm 

xo  ,.2  d2h(y)  ,m  q Exo  mq  Exo,  .dh(y)  [NDD(0)]  2 

2 ^ 2+'  22+  2 2'^  Hv  I y 

dy^  4kV  2kV  dy  KDn 


h(y) 


Rewrite  (2.8)  as 


y2dfhMt  . b2y2h(y)  = 0 


dy 


where  a = 1 + (2/m) , 


(2.5) 

(2.6) 

(2.7) 

= 0 . 

(2.8) 

(2.9) 

2.10) 


and  b2  = (2kT/raqExo)2[NDD(0)f/KDp 


(2.11) 
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and  b is  the  positive-valued  square  root  of  (2.11).  The  solution  of 
(2.9)  is  on  page  440  of  [20]: 

h(y)  > y"1/m  [Cjl^tby)  + C2K1/m(by)]  (2.12) 

where  and  C£  are  arbitrary  constants,  and  I^m  (K^)  is  the  1/m-th 
order  modified  Bessel  function  of  the  first  (second)  kind.  The  solution 
of  (2.6)  is  the  combination  of  (2.12)  and  (2.7): 

P(x)  = exp(-qExox/2kT){C1I1/m[bexp(mqExox/2kT)] 

+ C2Kl/m^bexp^mqExox/2kT^}  * (2.13) 

The  constants  and  C2  can  be  determined  by  the  boundary 
conditions  at  the  base  edge  of  the  space-charge  region  (x=0)  and  at  the 
back  surface  (x=W): 

P(0)  = [n2/NDD(0)]exp(qV/kT)  (2.14) 

5pExoP(W)  - DpC>  x=W  ■ sp(“)  • <2-15> 

In  (2.14),  as  in  (2.4)  previously,  we  have  assumed  the  full 
ionization  of  donor  atoms.  Using  (2.14)  and  (2.15)  with  m = 1 as  an 
example,  we  obtain 
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[BKj ( A)  + CK0(A)3[n?/NDD(0) ]exp(qV/kT) 

C1  = K'1(b)LBK1(A)  + CKQ(A)J  - K1(b)LBl1(A)  - CIQ(A)J  = C1  exP(qV/kT) 

(2.16) 


„ ^(A)  ~ CI0(A)][n2/NDD(0)]exp(qV/kT) 

2 = K^CblLBI^A)  - C1q(A)J  - I1(b)LBK1(A)  + CKq(A)J  = C2  exP(qV/kT) 

(2.17) 


where  A = bexp(qEx0W/2kT) , B = -S/0p , and  C = AqEx0/2kT. 

The  total  hole  current  density  is 

Jp  - Jp(0)  - q-ppExoP(0)  - qDp(f)x=0 

1/2  (2.18) 

= q[BpNDD(0)/K]  [CJIo(b)  - C'Ko(b)]exp(qV/kT)  , 

and  the  pre-exponential  factor  (often  called  the  reverse  saturation 
current)  is 

°P0  = <lC5pND0(0)/K31/2[c;io(b)  - CJK0(b)]  . (2.19) 

2.3  Exampl e 

The  foregoing  analysis  is  general,  subject  only  to  the  boundary 
conditions  (2.14)  and  (2.15).  We  now  consider  the  (retarding)  drift- 
field  (DF)  p/n  junction  Si  solar  cell  as  an  example.  Table  2.1  lists 
Jpo  for  different  values  of  W and  S.  Here  NDD(0)  and  NDD(W)  are  assumed 
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Table  2.1 

DF  Si  Solar  Cell  with  NDD(0)  = 1017  cm"3  and  NDD(W)  = 5 x 1018  cm'3 


W(pm) 

S(cms"l) 

Dp  (cm2s"1) 

JpQ  (10'^A/cm2) 

15 

106 

3.19 

6.82 

15 

103 

3.19 

2.27 

15 

0 

3.19 

1.75 

25 

106 

3.19 

5.31 

25 

103 

3.19 

3.15 

25 

0 

3.19 

2.77 

50 

106 

3.19 

5.23 

50 

103 

3.19 

4.77 

50 

0 

3.19 

4.64 

75 

106 

3.19 

5.90 

75 

103 

3.19 

5.80 

75 

0 

3.19 

5.76 

100 

106 

3.19 

6.54 

100 

103 

3.19 

6.51 

100 

0 

3.19 

6.51 

200 

106 

3.19 

8.10 

200 

103 

3.19 

8.10 

200 

0 

3.19 

8.10 
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to  be  10”  cm  3 and  5 x 10^  cm"3  respectively,  and  Dn  is  calculated  by 
W p 

using  Dp  = [ / Dp(x)dx]/W,  where 
o 

Dp(x)  = (47.7  + (495  - 47.7)/[l  + (N0[)(x)/6.3  x 1016) °‘76] }( kT/q)  [21]. 

For  comparison.  Table  2.2  lists  Jpo  for  a Si  solar  cell  with 
uni  form  base  doping  concentration  (1017  cm"3)  and  a back  ohmic  contact 
(BOC).  Note,  for  example,  that  JpQ  of  a DF  solar  cell  is  smaller  than 
that  of  a uniformly  doped  solar  cell  by  a factor  of  20  for  W = 25  pm. 
In  contrast  with  a DF  solar  cell,  a BOC  solar  cell  with  a thin  uniformly 
doped  base  region  (Lp  >>  W)  shows  a high  sensitivity  between  Jpg  and 
S.  Here  by  S we  mean  an  effective  surface  recombination  at  a low-high 
junction,  as  in  the  BSF  solar  cell.  The  DF  cell  has  the  advantage 

relative  to  the  BSF  cell  of  avoiding  patchiness  and  other  defects 
discussed  recently  by  Sah  and  co-workers  [22]. 

2.4  Discussion 

We  have  presented  general  analytic  results  for  JpQ  for 
nonilluminated  semiconductor  bipolar  devices  with  built-in  drift  fields 
under  the  assumption  that  minority-carrier  lifetime  has  a power-law 
dependence  on  the  shallow-level  donor  concentration.  The  assumption  is 
consistent  with  the  longest  lifetimes  measured  and  with  the  physical 
chemistry  of  di vacancy-donor  reactions  [15,18], 

As  an  example,  we  have  considered  a Si  DF  p/n  junction  solar  cell 
and  have  shown  that  (a)  Jpg  has  a nearly  constant  value  over  a broad 
range  of  W (from  15  to  100  pm)  if  the  back  surface  is  covered  by  ohmic 
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Table  2.2 

BOC  Si  Solar  Cell  with  Nqq  = 10^  cm-2  = constant 


W (pm) 

Dp  (cm2s_1) 

Jpg  (10"^A/cm2) 

15 

6.5 

14.6 

25 

6.5 

8.80 

50 

6.5 

4.57 

75 

6.5 

3.21 

100 

6.5 

2.58 

200 

6.5 

1.83 
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contact;  (b)  for  this  range  of  W,  Table  2.1  for  a DF  cell  shows  that 
Jpg  corresponds  to  large  open-circuit  voltage  as  calculated  from 
VqC  = V-pAn(Jgg/Jpg) : ^ ^or  ^SC  = and  667  mV  for  J<jg  = 
15  mA/cm2.  In  contrast,  from  Table  2.2  for  a BOC  cell,  for  comparison, 
VgC  = 590  mV.  Subject  to  the  considerations  discussed  below,  these 
values  of  Vgg  for  a DF  p+/n(x)  cell  with  W = 35  |jn  are  in  good  accord 
with  measured  VQC  (636  mV,  27°C)  in  [23],  the  first  VQC  reported  for  Si 
substantially  larger  than  600  mV,  and  with  the  Vgg  reported  for  a recent 
BOC  cell  (W=18  pm)  [24].  Note  that  the  lower  bound  of  this  range  of  W 
exceeds  the  average  light  penetration  depth  for  Si  for  the  AMI  solar 
spectrum.  Thus,  Jpg  from  the  quasi -neutral  base  imposes  no  significant 
limit  on  the  Vgg  achievable;  the  limit  will  be  established  by  the 
recombination  current  in  p/n  junction  space-charge  region  or  in  the 
quasi-neutral  emitter.  The  thickness  W,  however,  will  influence  Jgg. 
Because  increasing  W increases  the  internal  quantum  efficiency  while 
decreasing  Ex0,  we  may  expect  that  the  maximum  achievable  EFF  will  be 
obtained  over  a broad  range  of  W.  This  implies  that  a thin  base  of  a DF 
cell  can  save  material  and  fabrication  cost  at  no  loss  in  EFF. 

Two  approximations  used  in  his  work  deserve  further  emphasis.  We 
have  assumed  a minority-carrier  lifetime  that  varies  in  accord  with 

(2.1) ,  in  which  K,  for  m = 1,  is  given  numerically  directly  below 

(2.1) .  This  approximation  corresponds  to  the  longest  lifetimes  observed 
[18].  Thus  the  results  in  our  example  of  Sec.  2.3  correspond  to  best- 
case  solar-cell  performance  for  the  drift-field  structure  studied.  As  a 
second  approximation,  we  have  assumed  an  exponential  dependence  on  x of 
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the  donor  density;  for  full  ionization  of  the  donor  atoms,  this  gives  a 
position-independent  built-in  electric  field.  In  practice,  one  does  not 
see  exponential  profiles  of  donor  atoms,  but,  rather,  in  general,  a more 
complicated  mathematical  function  (not  the  Gaussian  or  erfc  functions 
commonly  assumed  either).  The  equilibrium  electric  field  becomes 
further  complicated  by  the  compensation  arising  from  the  p/n  junction 
formation.  The  sensitivity  to  doping  profile  of  the  results  of  Sec.  2.3 
deserves  further  study;  but  spot  checks  using  the  numerical  methods  and 
profiles  of  [25]  suggest  the  good  accuracy  of  our  results. 


CHAPTER  3 

MEASUREMENTS  AND  MEASUREMENT  ANALYSIS 
OF  NON-OPTIMAL  SOLAR  CELLS 


3.1  Introduction 

In  Chapter  2,  we  have  shown  theoretically  that  a solar  cell  having 
a drift-field  base  region  can  have  low  dark  current  and  thus  high  open- 
circuit  voltage.  We  now  study  experimentally  the  performance  of  DF 
solar  cells. 

We  investigate  p/n/n(x)/n/n-substrate  solar  cells  made  for  this 
study.  These  cells  have  an  n/n(x)/n  base  region  in  which  the 

n(x)-dri ft-fiel d region  has  a thickness  of  11  pm.  Comparison  with  the 
analysis  of  Chapter  2 shows  that  these  cells  are  highly  non-optimal . 
But  the  analysis  to  follow  on  the  non-optimal  cells  has  general 
applicability,  and  thus  is  intended  to  benefit  studies,  planned  for  the 
future,  of  DF  cells  having  thicker  retardi ng-fiel d regions. 

We  aim  here  to  estimate  minority-carrier  lifetime  in  the  base 
region  first  from  the  measured  short-circuit  current.  Second,  we  shall 
use  the  method  of  ionizing  radiation  due  to  Rosenzweig  [26],  Because 
the  non-optimal  devices  exhibit  performance  consistent  with  the 
existence  of  sizable  shunt  conductance,  we  employ  measurements  of  the 
short-circuited  device  to  avoid  the  influence  of  this  conductance. 
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We  also  qualitatively  determine  the  predominant  components 
contributing  the  dark  current.  The  qualitative  determination  comes  from 
employing  the  non-il  1 uminated  I-V  characteristics  together  with  the 
measured  fill  factor  and  spectral  response.  Use  of  the  estimated  base 
lifetime  yields  a semi-quantitative  determination. 

Finally,  we  explore  implications  deriving  from  comparing  the 
performance  of  the  non-optima!  DF  cells  to  that  of  extended-BSF  cells 
analyzed  by  computer  [27]. 


3.2  Device  Make-Up 

The  cells  used  here  were  fabricated  on  antimony-doped  substrates. 
Substrate  resistivity  is  approximately  0.01  ohm-cm  and  orientation  is 
three  degrees  off  the  (111)  direction.  Epitaxial  layers  of  Si : P 
chemically  vapor  deposited  (CVD)  on  the  substrates  constituted  the 
n/n(x)/n  base  region.  The  n-region  nearest  the  surface  is  approximately 
uniformly  doped  and  has  a doping  density  of  6 x 1016  cm"3  and  a 
thickness  of  9 pm.  The  doping  density  in  the  n(x)-dri ft-f ield  region 
(thickness  11  pm)  increases  exponentially,  to  a good  approximation,  from 
6 x 1016  cm  3 to  3 x 10^3  cm"3,  which  is  the  doping  density  of  the 
uniformly  doped  n-region  (thickness  18  pm)  nearest  the  substrate 
(Fig.  3.1).  After  the  epitaxial  growth,  a shallow  p+  emitter  (0.2  pm) 
resulted  from  boron  diffusion.  A metal-grid  pattern  providing  8% 
shading  was  used  together  with  a mul  ti -1 ayer-anti refl ection  (MLAR) 
coating.  All  cells  were  cut  to  form  a square  of  approximately  2 cm  x 


2 cm. 


CARRIER  CONCENTRATION  (cm- 


a kj 


Fig.  3.1 


The  doping  density  for  epitaxial  region 
obtained  by  spreading  resistance  measurement 
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3.3  Measured  Cell  Performance  and 
Qualitative  Interpretation 

The  performance  of  the  DF  solar  cells  was  measured  for  an  incident 
AMI  solar  spectrum  with  a radiant  power  density  of  100  mW/cm2  at  28°C. 
Table  3.1  shows  fill  factor  (FF),  short-circuit  current  density,  open- 
circuit  voltage,  and  conversion  efficiency  for  four  such  cells. 

Figure  3.2  is  the  dark  I-V  characteristic  for  the  drift-field 
cell.  The  current  is  proportional  to  exp(qV/mkT)  having  m=2  for  a 

large  range  of  V but  nowhere  is  m < 2.  Thus  graphical  methods  [28]  for 
separating  the  quasi-neutral -region  current  from  the  space-charge-region 
current  cannot  be  employed. 

To  more  closely  examine  the  quality  of  the  junction,  we  measured 
the  I-V  characteristics  of  many  small -area  diodes  made  from  the 
2 cm  x 2 cm  solar  cell.  To  make  such  cells,  we  cover  the  entire  front 

surface  with  aluminum,  putting  wax  on  the  back  surface  and  some 

randomly-chosen  spots  of  the  front  surface.  Then  we  etched  away  the 
aluminum  layer,  the  contact  fingers,  the  anti  reflection  layer  and 
through  the  p+  front  layer  in  the  area  not  protected  by  wax.  This 
procedure  yields  small  diodes  whose  front  contacts  are  isolated  from  one 
another. 

We  measured  I-V  characteristics  of  eleven  such  diodes  whose  area 
ranged  from  1.82  x 10-2  cm2  to  8.3  x 10"2  cm2.  Figure  3.3  shows  two  of 
them.  Curve  (1)  behaves  more  like  a normal  diode  (m  = 1 for  some  range 
of  V)  while  curve  (2)  deviates  far  from  normal  diode  behavior  (m  > 2 for 
all  V).  The  other  nine  curves  fall  between  these  two  extreme  cases. 
Thus,  the  areal  inhomogeneity  of  the  junction  contributes  to  its  poor 

qual  ity. 
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Table  3.1 


Parameters  of  the  drift-field  solar  cells 
under  AMI,  28°C,  100  mW/cm2 


V0C 

JSC 

FF 

EFF 

(mV) 

(mA/cm2) 

(%) 

4360-2-1 

596 

24.5 

0.77 

11.2 

4360-2-2 

596 

24.5 

0.79 

11.5 

4360-8-1 

590 

24.7 

0.67 

9.7 

4360-8-2 

598 

25.0 

0.78 

11.6 

I (A) 
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Fig.  3.2 


Dark  I-V  characteristic  for  the  drift-field  cell. 

The  current  is  proportional  to  exp(qV/mkT) , where  m > 2. 


I (A) 
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Fig.  3.3 


Dark  I-V  Characteristics  for  two  small -area  diodes 
made  from  a 2 cm  x 2 cm  drift-field  cell. 


This  areal  inhomogeneity  can  come  from  several  different  origins. 
By  areal  inhomogeneity  we  mean  any  variation  in  the  junction  properties 
occurring  across  the  area  of  the  solar  cell,  which  is  a large-area 
device.  Examples  include  propagation  of  structural  defects  on  the 
surface  into  the  device  during  thermal  processing,  which  can  lead  to 
localized  portions  of  the  area  in  which  shunt  conductance  is  high,  as 
pointed  out  first  by  Stirn  [29];  leakage  at  the  rim  of  the  solar  cell, 
again  giving  rise  to  high  shunt  conductance;  inhomogeneities  in  the 
distribution  of  impurity  centers,  either  in  the  junction  space-charge 
region  or  in  the  quasi -neutral  regions  [30,  31];  inhomogeneity  across 
the  front  surface  of  the  solar  cell  arising  from  the  mixture  of  ohmic- 
contact  grid  lines  with  AR  coating  or  thermal  oxide  [30,  31]. 

The  good  fill  factor  listed  in  Table  3.1,  also  shown  in  Fig.  3.4, 
rules  out  the  possibility  of  the  dominance  of  p/n  junction  space-charge- 
region  recombination  current  throughout  the  entire  volume  of  the  space- 
charge  region.  It  does  not  rule  out,  however,  an  areal  inhomogeneity  in 
this  space-charge  region  current,  and  the  presence  of  such  inhomogeneity 

is  consistent  with  results  previously  discussed  and  illustrated  in  Fig. 

3.3. 

In  general  we  may  expect  the  value  of  shunt  conductance  to  depend 
on  V,  particularly  if  surface  channels  on  the  rim  of  the  solar  cell 
contribute  to  this  conductance.  Nonetheless,  consistently  with  the 
conclusions  of  several  authors  [32,  33],  we  may  qualitatively  expect 
that  shunt  conductance  can  result  in  a reciprocal  slope  m > 2.  From  the 


I (mA) 
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V (VOLTS) 


Fig.  3.4 


Illuminated  I-V  characteristics  for  drift-field 
cell  obtained  under  AMI,  28°C,  100  mW  crrf2. 


reverse-bias  I(V)  characteristics,  we  estimate  that  shunt  conductance 
Ggiy-  2.5  x 10“^  mhos/cm^. 

We  now  consider  implications  of  the  spectral  response  in 
Fig.  3.5.  The  DF  cells  suffers  from  poor  red  (X  > 0.9  pm)  response  and 
poor  blue  (x  < 0.5  pm)  response.  The  poor  red  response  might  result 
from  either  the  thinness  of  the  base  region  in  which  the  retarding  drift 
field  exists  or  from  low  minority-carrier  lifetime  in  the  base  region. 
The  poor  blue  response  might  be  caused  by  low  minority-carrier  lifetime 
in  the  emitter  region  or  high  recombination  velocity  at  the  front 
surface  or  both. 

The  assorted  foregoing  evidence  suggests  that  the  dark  current 
coming  from  the  areal  inhomogeneity  discussed  above  and  from  the  quasi - 
neutral  emitter  and  base  regions  may  be  important  in  limiting  the 
performance.  But,  as  yet  we  have  advanced  no  evidence  to  suggest  which 
quasi -neutral  component  of  current  dominates. 

For  closer  examination  of  the  influence  of  each  component  of  the 
dark  current,  we  turn  to  the  following  estimate  of  the  mi nority-carrier 
lifetime  in  the  base  region  in  a one-dimensional  model. 

3.4  Estimation  of  Minority-Carrier  Diffusion 
Length  and  Lifetime  in  the  Base  Region 

We  now  try  to  estimate,  by  theory  and  experiment,  the  minority- 
carrier  diffusion  length  and  lifetime  in  the  base  region.  First,  using 
the  simplifying  assumption  of  position-independent  lifetime  we  derive  an 
expression  for  short-circuit  current  density  in  which  the  diffusion 
length  Lp  appears  as  a parameter.  This  expression  yields  a family  of 


so 
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(m/v)  asNOdsaa  ivaxoads 


Spectral  response  for  the  drift-field  cell. 


27 


curves  in  this  parameter;  comparison  with  the  measured  Jsc  then 
determines  Lp.  Second,  we  use  the  method  of  X-ray-induced-photocurrent 
to  estimate  the  diffusion  length  and  lifetime. 

3.4.1  Position-Independent  Lifetime  and  Measured  J^p 

Let  the  quasi -neutral  portion  of  the  uniformly-doped  n-region 
nearest  the  surface  region  in  Fig.  3.1  be  called  region  1 and  have 
thickness  W^.  Let  the  n(x)-drift-field  region  (as  region  2)  have 
thickness  W2  and  the  region  nearest  the  substrate  of  x-i ndependent 
doping  concentration  (region  3)  have  thickness  W3.  Let  the  thickness  of 
the  quasi -neutral  p+  emitter  together  with  the  space-charge-region 
thickness  be  WQ. 


For  AMI  illumination  applied 

to 

this  solar  cell. 

the  following 

equations  hold: 

.2 

d Pi  p, 

Dp  h2  - * * -SE<X> 

K d x p 

i n 

region  1 

(3.1) 

d2p2  dp?  p„ 

i n 

region  2 

(3.2) 

.2 

d p3  p3 

°p  A - i - 

in 

region  3. 

(3.3) 

Here  G = [An(N3/N1)]/(W2  - W1), 

(3.4) 

in  which  Nj  and  N3  are  the  uniform  phosphorous  doping  concentrations  of 
regions  1 and  3 respectively. 
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The  boundary  conditions  are 


Pl(W0)  = 0 (3.5) 
Pl(Wj)  = P2 ( ^1 ) (3.6) 
Jl(Wj)  - J2(Wj)  (3.7) 
P2(W2)  = P3(w2^  (3.8) 
J2(W2)  = ^3(^2^  (3.9) 
J3(W3)  = qSp3(W3)  . (3.10) 


Here  S is  the  surface  recombination  velocity  at  the  interface  between 
the  epitaxial  region  and  the  substrate. 

To  model  the  optical  hole-electron  generation  rate  g^(x)  as  a 
function  of  x,  we  use  the  approximation  of  Hsieh,  Hu,  and  Drowley  [34]: 

4 g . 

9F(x)  = l x A • (3.11) 

t j=1  x+  aj 

Here  g^,  g3,  g3,  g^,  a^,  a2»  a3,  and  a^  are  3.7  x lO1^,  77  x 1013, 
-77  x 1015,  4.95  x 1015  cm'2  s'1,  6.37  x 10"7,  5.6  x 10'5,  6.14  x 10"4’ 
and  2.23  x 10~3  cm  respectively.  The  solutions  for  (3.1),  (3.2)  and 

(3.3)  are 
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y v L 4 x+a. 

P1  = AsexP(— ] + Bsexp(-  — } + ^ l qiK{-TA)  > 

P P P J=1  p 


p2  = Cs 


, x,  x (Ln  Lp)  4 X+a. 

exp(^)  + D exp(-  — ) + 2- j 9n-K(+-Ti) 

Ln  L-  D„(L>L-)  j=l  J t* 


P'  P P' 


and 


L„  4 x+a. 


p3  Esexp(|_  ^ + Fsexp(“  + TF-  j gjK(  L ^ ’ 


P J=1 


respectively.  Here  As,  B$,  Cs,  Ds,  E$  and  F,.  are  constants,  and 


_1  = [(i) 2 + M/2  ^ 

Lt  u2;  ^2J  +2  , 


x+  a . x+  a . x+  a . x+  a . x+  a . 

Kf-i-J-)  = exp(-  -rJ-)Ei(-rA)-exp(-rJ-)E1(-  -j-J-)  , 

P P P P p 


and 


x+a,  x+a.  x+a.  x+a.  x+a. 

K("— T-)  = exP( -)M -)-exp( — r^-)E,  ( -J-) 

lZ  L 1 L"  L 1 L 


in  which  Ei  is  the  exponential  integral  function.  The  short 
current  density  is 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


ci rcuit 


JSC  " "qDp(Hf)x=Wf 


x+a . , x+a . 

Here  H(-j-l)  = K(-j-J-) 

P P 


x+a . 


x+a . 


1 /\'U.  A'  a . X + 3.  A Ta. 

— [exp(-  -rJ-)Ei(-^)  + exp(-j-J-)E.(-  -j— J-) 


x+a . 


Boundary  conditions  (3.5)  through  (3.10)  determine  constants  A 


and 


Lo  Wn  4 wn+ai  2wn 

Bs  “ 20^  exP(— ) l giK(  V)  "As  exP(-r}  * 
P P J=1  J 0 S D 


P P J 


Here 


sn 


sA 

, W?  5T~  2W,-W, 

R3  + {]H!exp(— ) + — ^exp( — r- — -)  ] 

P P S-^  P 


(3.18) 


(3.19) 

s and  Bs: 

(3.20) 

(3.21) 
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D 

S-K-O. 

W~  3 L 2W--W,  , 

■[exp(-r^) AxpC-g  2)](G  + -£)} 

P S->  P Lo 


W2+WL 


x exp(— p-)  [-(G  - — + r^) 
LP  LP  LP 


L Wn-W.  4 W +a . 

2D^XP(  -L--~)  I 9iK(~LL)  + RJ 
P P j=l  J P 


D 

Sh — 2. 

,1  W9  Ln  2W--W, 

- {f-[exp(-p)  +— p£exp(-[_  2)] 

p p s__e.  p 
lp 


■[exp(-p) 

P 


D 

S 

Ld  2Vw?  i -Vwi 

P-  exp(  ^ -)  ](G-  -*■)}  exp( — - — - 

l P LP  LP 


S-  ^ 
LP 
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and 


sd 


0 

c,  P 

, W?  aT“  2W,-W, 
-{f" [exp(-j-)  + — ^ exp(~L- -2-)  ] 

p P s-J-  P 


D 

s+-£ 

W9  3 L 2W~-W9  - 

■[exp(-j-=-) ^ exp(  -2')](G  + -}•)} 

P SA  LP  Lp 


VW1  i i wi  i i 2Wn-W1 

x exp(— =— ^) [ (G  - 4 - T^expM-MG  - -±-  + 4)exp(-f-4 ] 
k k Lp  Lp  Li  Lp  Lp 


■p  “p  P P ~p 


D 

— 2. 

i w,  2W~-W9 

+ {f— [exp(— ) + — p2-  exp( — r~h] 

P P S-^  P 


D 

C_L P 

w9  L 2W--W, 

[exp(-j^) 4 exp(— P---)](G 
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4)1 

LP 
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LP  LP  LP 
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P Lp  p 
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The  following  three  constants  in  (3.22)  are  given  by: 


4 g,L 


W,  + a . W.  + a . 


Ri  ■ i w - - h(V)] 

J=1  P p Lp 


4 gf(L*  L”) 


W.+a.  W.+a. 


, ,3,4, 


P P P 


4 9,L  i W +a . W.+a. 

R2  = X ~2D^  [<G  +4)K(Jrl)  - Hf-t-J-)] 
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+ [exp(^) 
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Here  H(;— ^L)  = J*K + GK^-J-) 
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= (G j)  exp( J-)Ei( *4 

LP  LP  LP 


x+a . 


1 A 1 a . x+  a . 

- (G  + “■)  exp( I^E.  ( T1) 

Ln  L 1 L 

p p p 


(3.27) 


Table  3.2  lists  J$c,  which  is  insensitive  to  S,  for  different 
values  of  diffusion  length  Lp.  Comparing  Jsc  in  Table  3.2  with  the 
measured  Jsc  in  Table  3.1,  yields,  in  the  present  model,  the  minority- 
carrier  diffusion  length  Lp  is  about  20  pm.  The  diffusivity  Dp  is 

O 

assumed  to  be  4.0  cnr/s,  consistently  with  the  data  reviewed  in  Muller 
and  Kami  ns  [21]  for  a phosphorus  concentration  of  about  3 x 1017  cm"3. 
Thus,  the  minority-carrier  lifetime  is  about  1 ps  obtained  from 

tp  = lp/dp- 
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Table  3.2 

The  theoretical  variation  of  short-ci  rcuit  density  Jor.  with  minority- 
carrier  diffusion  length  Lp  under  the  assumption  of  position  independent 
lifetime  in  the  non-uni formly  doped  base  region  of  a DF  solar  cell. 


LP 

(pm) 

J$c  (mA/cm^) 
S = 10  cm/s 

S = 10®  cm/s 

1 

7.6 

7.6 

5 

16.3 

16.3 

10 

21.2 

21.2 

20 

25.1 

24.6 

30 

26.3 

25.5 

40 

26  .8 

25.9 

50 

27.0 

26.0 

60 

27.2 

26.1 
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If  the  parameters  have  the  values  determined  by  this  first  method, 
then  the  lifetime  is  long  enough  that  a similar  solar  cell,  whose  entire 
epitaxial  base  region  is  a drift-field  region  of  about  50  microns,  will 
have  a conversion  efficiency  exceeding  15%,  according  to  the  exact 
computer  solutions  reviewed  by  Lindholm,  Neugroschel  , and  Sah  [35], 
Thus  it  becomes  of  interest  to  seek  other  evidence  to  support  the 
results  obtained  by  this  method.  We  now  turn  to  a method  of  measurement 
suggested  by  Rosenzwig  [26]. 

3.4.2  X-ray  Induced  Photocurrent 

We  expose  the  drift-field  solar  cell  and  a standard  cell  with  known 
diffusion  length  to  low  energy  X-ray  irradiation  and  measure  the  X-ray- 
generated  of  these  two  cells  simultaneously.  It  has  been  shown  [36] 
that  the  carrier  generation  rate  of  X-ray  irradiation  is  quite  uniform 
and  can  be  treated  as  spatially  constant.  The  constant  generation  rate 
gE  can  be  obtained  by  substituting  the  measured  Jsc(2  x 10"8Acm"2) , the 
diffusion  length  (Lp  = 157  urn),  and  the  base  width  (Wg  = 320  urn)  of  the 
standard  cell  into  the  following  equation, 

WR 

coshj-2-  - 1 

Jsc  = 19ELp  Ejj . (3.28) 

sinh^ 

P 

The  gE  = 1013  cm'3. 

Because  the  base  region  of  the  standard  cell  is  uniformly  doped, 
(3.28)  can  be  obtained  by  solving  (3.1)  with  boundary  conditions 
P^(0)  = 0 at  the  base  edge  of  the  space  charge  region  and  P^(Wg)  = 0 at 
the  back  surface. 
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To  find  the  short-circuit  current  density  for  the  DF  cell  for 
constant  generation  rate,  we  solve  (3.1)  through  (3.10)  and  now  use 
constant  generation  rate  g^  instead  of  the  position-dependent  generation 
rate,  (3.11).  Solutions  (3.18)  through  (3.26)  can  be  used  after  some 


minor  modifications: 


L 

2D 


Vii 


L 4 W +a . 

delete  2^-  z g.H(-p — i)  in  (3.18);  change 

P j=l  J P 


L ) 1n  ^.21)  an^  (3.22)  to  g_L  /D  ; rewrite  (3.24), 

j=l  J p t p p 

1 


(3.25)  and  (3.26)  as 


R1  ■ G9eVD( 


(3.29) 


R2  * G9e7DP 


(3.30) 
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^£exP(-f^)][^Li 

S—P.  P P 
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^exp( 
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(3.31) 
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Substituting  measure  Jsc  (1.6  x 10'8  A cm'2)  for  the  DF  cell  into 
modified  (3.18)  through  (3.26)  yields  that  Lp  is  about  16  pm  and  Tp  is 
about  0.65  ps. 

The  results  are  close  to  those  obtained  by  the  estimation  of  the 
last  section:  Lp  is  about  20  pm  and  tp  is  about  1 ps. 

3.5  Components  of  the  Total  Non-Il 1 umi nated  Current 
The  treatment  of  the  foregoing  section  has  determined  the  minority 
carrier  lifetime  t for  the  quasi-neutral  base  region  of  the  non-optimal 
drift-field  solar  cell  in  a one-dimensional  model  that  has  assumed 
position  independence  of  t.  We  now  use  this  result  to  calculate  the 
reverse  saturation  current  (pre-exponential  factor)  of  the  quasi -neutral 
base.  By  exploring  its  relation  to  the  measured  open-circuit  voltage 
and  short-circuit  current,  we  can  then  estimate  the  additional  current 
components  that  contribute  to  the  values  of  these  observables.  All  of 
the  quantitative  treatment  in  this  section  assumes  a one-dimensional 
model;  we  have  already  discussed  the  inadequacy  of  such  a model  because 
of  the  presence  of  various  mechanisms  that  contribute  to  areal 
inhomogeneity  and  lead  to  a reciprocal  slope  m > 2 in  the  log  (I)  vs. 
qV/kT  measured  characteristic . Thus  the  origin  within  the  device  of  the 
additional  current  components  referred  to  above  cannot  be  precisely 
determined. 

The  continuity  equations  for  the  three  portions  of  the  quasi - 
neutral  base  region  of  the  non-optimal  DF  solar  cell  are  as  follows: 
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d P1  P1 
D — i - — = 0 

P dx  Tp 


in  region  1 


(3.32) 


d2P-  dP9  P9 

D i + G -4  - -L  = 0 

p dx2  dx  Tp 


in  region  2 


(3.33) 


d P3  P3 

Dn  i - -i  = 0 

P dx  Tp 


in  region  3 


(3.34) 


The  six  boundary  conditions  are 


together  with  (3.6)  through  (3.10).  The  solutions  for  (3.32)  - (3.34) 
are 


= Aexp(x/Lp)  + Bexp(-x/Lp) 


(3.36) 


P2  = Cexp(x/Lp)  + Dexp(-x/L~) 


(3.37) 


and 


P3  = Eexp(x/L  ) + Fexp(-x/L  ) 

■j  p p 


(3.38) 


respectively. 
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To  determine  Jpo, 


= "<1(l£’)(a  - B)/exp(^) 


(3.39) 


we  need  only  to  determine  constants  A and  B: 


A " exp^kT^  a" 


(3.40) 


and 


B " exp(kT)  " A 


(3.41) 


Here 
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, ir  "o  2VL-VL  H 

{Hexp(— ) + exp(  ^ -) ] - (G  +-l)[exp(r^)  -exp(— 2-1)]} 

p P P L_  P ld 


W„  2W„-Wr 


W?  , . w,  W, 

x exp(— )(G  — - + — )exp( )exp( — j) 

LP  LP  LP  LP  Lp 


r lr  2W--W?  W,  2W~-W0 

-{[H^exp^)  + exp(  - - -)]-  (G  - ^)[exp(^)  - exp(  - 2)]} 
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W,  W, 
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(3.42) 
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(3.43) 


From  the  two  methods  discussed  in  Sec.  3.4,  we  determined  the  values  of 
the  lifetime  and  the  diffusion  lengths  to  be 


Tp  = 1 ns,  Lp  = 20  pm  (method  1) 


and 


Tp  = ^*65  Ps»  ~Lp  = 16  pm  (method  2) 


For  the  reverse  saturation  current,  these  values  yield 


Jn_  = 6.0  x 10 
po 


-13 


Acm 


-2 


(method  1) 


and 


Jpo  = 9.0  x 10 


-13 


Acm 


-2 


(method  2) 


Use  of  the  idealized  expression  [37]  for  the  open-circuit  voltage. 


VQC  = 625  mV  (method  2) 

This  contrasts  with  the  measured  open-circuit  voltage  of  598  mV  for  the 
cell  under  study. 

These  discrepancies  suggest  the  following  interpretations.  One 
possibility  is  that  the  areal  inhomogeneity,  perhaps  particularly  the 
shunt  conductance,  is  responsible  for  the  discrepancy.  The  high  fill 
factor  (FF  = 0.78)  rules  out  the  junction  space-charge-region 


(3.44) 


we  get 


VQC  = 636  (method  1) 


and 
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recombination  current  as  a dominant  contributor  in  a one-dimensional 
model.  Another  possibility  is  that  the  recombination  current  in  the 
quasi-neutral  emitter  also  contributes  to  the  discrepancy;  the  poor  blue 
response  shown  in  Fig.  3.4  supports  this  conclusion. 

3.6  Further  Comparison  with  Theory 
The  theoretical  results  of  Chapter  2 were  based  on  the  assumption 
of  a best  lifetime  in  Si  for  a given  doping  concentration.  Thus 
T = f(x)  f°r  the  Dfr  cell  studied  in  that  chapter.  A similar  study,  by 
Sah  and  Lindholm  [27],  assumed  a best  lifetime  for  Si  for  an  extended 
back-surface-field  in  which  the  drift-field  region  penetrated  various 
distances  into  the  n-type  base;  large  penetrations  corresponded  nearly 
to  a drift-field  solar  cell  and  small  penetrations  corresponded  to  a BSF 
cell,  thus  enabling  a systematic  comparison  between  the  two  types  of 
cells.  The  Sah-Lindholm  work  solved  the  finite-difference  counterparts 
of  the  seven  Shockley  equations  by  a contact -to-contact  numerical 
solution  for  a cell  thickness  of  50  pm. 

The  results  reported  in  Chapter  2 agree  closely,  for  this 
thickness,  with  those  obtained  by  Sah  and  Lindholm  for  cells  having 
large  BSF  penetration  for  an  assumed  concentration  of  Zn  atoms 
(recombination  centers)  of  1012  cm'3,  which  corresponds  to  a best 
lifetime  in  n-type  silicon  for  a donor  dopant  concentration  of 
^ cm  • This  dopant  concentration  was  assumed  as  the  lowest  base 

dopant  concentration  present  in  the  extended  back-surface-field  cell 
treated  by  Sah  and  Lindholm.  This  Zn  concentration  corresponds  to  a 
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minority-carrier  lifetime  of  about  20  ps,  which  agrees  with  the  best 
lifetime  assumed  in  Chapter  2. 

The  Sah-Lindholm  calculations  predict  an  open-circuit  voltage  of 
640  mV,  for  x = 1 ps,  and  630  mV,  for  x = 0.65  ps.  These  results  apply 
for  a cell  having  a total  thickness  of  50  pm  and  a back-surface-field 
penetration  of  40  pm,  which  is  the  structure  treated  by  Sah  and  Lindholm 
that  is  closest  to  the  experimental  structure  we  have  studied  in  this 
chapter.  Further  Sah  and  Lindholm  assumed  no  metal  shading  and  no 
reflection  of  incident  light  from  the  front  surface. 

If  we  consider  the  structural  differences  between  the  modeled  and 
the  experimental  solar-cell  structures,  we  see  that  the  conclusions  of 
the  previous  section  (Sec.  3.5)  still  apply:  most  probably  areal 
inhomogeneity,  including  the  presence  of  shunt  conductance,  together 
with  some  possible  degrading  influence  from  the  highly  doped  emitter 
region  explain  the  discrepancy  between  the  open-circuit  voltage 
calculated  from  a one-dimensional  model  of  the  quasi -neutral  base  region 
and  that  observed  experimentally. 

The  design  implications  are  clear.  Improved  epitaxial  growth,  both 
in  doping  profile  and  in  general  quality,  are  needed  to  improve  the 
quality  of  the  DF  solar  cell.  That  such  improvements  are  possible 
follow  from  the  work  of  D'Ailleo,  Robinson  and  Kressel , who  reported  in 
1975,  an  open-circuit  voltage  at  29°  C of  636  mV  and  a conversion 
efficiency  (corrected  for  excessive  metal  shading  and  reflection)  of 
about  15%  for  a thickness  of  35  pm  [23,38]. 
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3.7  Recent  Experimental  Findings 

Lindholm  and  Sah  [39]  have  reported  and  interpreted  experimental 
results  for  DF  cells  of  40  ^un  having  an  improved  doping  profile  relative 
to  that  described  in  the  earlier  sections  of  this  chapter.  The  measured 
performance  parameters  were  slightly  better  than  those  reported  above. 
A key  finding  is  that  the  drift-field  cells  showed  considerably  better 
performance  than  control  samples  made  from  thick  (12  mil)  Cz  silicon. 

The  control  samples  normally  would  exhibit  much  better  performance 
than  that  measured  for  the  OF  cells.  This  suggests  an  absence  of  good 
control  in  the  formation  of  the  Si :8  emitter  by  solid-state  diffusion  of 
boron,  because  this  formation  involves  the  only  processing  common  to  the 
DF  cells  and  the  control  cells. 


CHAPTER  4 

A HEURISTIC  APPROACH  TO  STEADY-STATE  SEMICONDUCTOR 
FILAMENT  THEORY,  WITH  APPLICATIONS 

4.1  Introduction 

The  broad  aim  of  this  study  is  to  derive  approximate  solutions  for 
time-independent  boundary-value  problems  that  arise  in  many  different 
fields  of  physics.  The  problems  feature  injection  of  a variable  of 
interest  at  one  planar  surface  combined  with  mechanisms  (such  as 
recombination)  that  tend  to  produce  a vanishing  of  this  variable  within 
the  volume  and  at  the  other  bounding  surfaces  including  the  lateral 
surfaces  normal  to  the  plane  of  the  injecting  surface.  We  illustrate  a 
method  of  solution  for  such  problems  by  specific  semiconductor 
examples.  One  such  example  appears  in  Sec.  4.2;  others  are  in 

Sec.  4.4.  In  these  examples  the  minority  carriers,  which  are  either 
holes  or  electrons,  are  injected  into  the  semiconductor  over  a barrier 
in  response  to  an  applied  potential  difference  V.  Thus  a characteristic 
exp(qV/kT)  source  term  arises.  Following  the  injection,  the  minority 
carriers  flow  by  diffusion,  and  recombine  in  the  volume  and  at  the  back 
and  lateral  surfaces,  where  the  surface  recombination  velocity  S is 
taken  to  be  infinite.  Thus  the  solution  applies,  for  example,  to 
semiconductor  p/n  junction  or  Schottky-barri er  devices  containing  grain 
boundaries,  with  infinite  S on  the  grain  boundaries,  as  in  [7], 
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Although  the  time-dependent  (transient)  counterpart  of  these  three- 
dimensional  boundary-value  problems  has  received  considerable  attention 
in  the  past,  the  time-independent  problem  of  determining  the  current 
(or,  in  general,  flux)  response,  I(V),  has  been  neglected.  As  we 
demonstrate,  the  exact  solution  of  this  problem  runs  into  mathematical 
difficulty.  Thus  much  attention  is  devoted  to  a method  of  approximate 
solution,  which  is  a heuristic  method  in  the  sense  that  it  provides  aid 
and  direction  toward  solving  a problem  though  full  rigorous 
justification  is  absent. 

Shockley  first  studied  this  class  of  semiconductor  problems,  which 
has  come  to  be  known  as  filament  theory  [6].  Here  we  limit 
consideration  to  quasi-neutral  filaments,  having  position-independent 
doping  concentrations,  that  obey  the  low-injection  condition.  Thus  the 
governing  differential  equation  is  the  diffusion  equation. 
Recombination  at  the  lateral  surfaces  produces  a three-dimensional 
boundary-value  problem;  volume  recombination  is  also  considered.  In 
heat  conduction  [40]  and  in  electromagnetism  [41],  analogous  problems 
arise,  although  the  analog  to  volume  recombination  vanishes.  Solutions 
for  this  class  of  problems  involve  linear  combinations  of  infinitely 
many  eigenfunctions. 

Time-dependent  solutions  for  excess  carriers  in  a semiconductor 
filament  have  been  found  by  many  researchers  [6,  42-44].  These  previous 
workers  have  determined  that  truncation  of  the  infinite  series,  keeping 
only  the  first  term,  serves  as  an  adequate  approximation. 
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When  this  approximation  is  applied  to  the  corresponding  time- 
independent  (steady-state)  boundary-value  problem,  one  can  write 
expressions  for  current  density  and  total  current,  which  are  variables 
of  interest  in  semiconductor  devices.  But  the  exact  solution  for 
current  gives  rise  to  mathematical  difficulties,  and  the  quality  of  the 
first-term  approximation,  used  for  the  transient  counterpart  of  these 
problems,  becomes  questionable. 

In  this  study,  the  origin  of  the  mathematical  difficulties  is 
described  in  a related  and  simpler  two-dimensional  boundary-value 
problem  for  which  a closed-form  solution  suffices.  We  modify  the 
boundary  conditions  to  eliminate  such  difficulties  and  specify  the 
extensions  needed  to  accomodate  the  three-dimensional  problem  of 
interest.  Thus  we  obtain  an  analytic  expression  for  the  current.  In 
this  expression  the  filament  lifetime  appears.  The  development  leads 
to  the  definition  of  this  parameter  in  both  first  and  higher  approxima- 
tions (that  is,  higher  modes  are  included  in  the  terminology  of  [41]  and 
[42]. 

4.2  Analysis 

4.2.1  First-Term  Approximation  for  Transients 

The  lifetime  of  excess  carriers  in  a filament  can  be  defined  in 
terms  of  the  decay  of  photoconductivity  in  a sample  of  known  geometry. 
Consider  a sample  in  the  form  of  a rectangular  parallel  piped,  bounded  by 
the  planes  x = 0,  x = 2A;  y = 0,  y = 2B;  z = 0,  z = 2C.  It  is  assumed 
that  the  sample  is  p-type  material  and  that  there  is  negligible  electric 
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field  acting  on  the  minority  electrons  in  the  sample,  which  is 
consistent  with  a filament  in  low  injection  having  a position- 
independent  doping  concentration.  Suppose  the  filament  has  been 
illuminated  for  t < 0;  the  differential  equation  governing  the  excess 
carrier  density  n is 


M = 0 = gE  + 0 v2n  , t < 0 (4.1) 
n 

where  g^  is  the  resulting  position-independent  generation  rate 
throughout  the  sample,  x is  the  minority  electron  lifetime,  and  D is 
the  electron  diffusivity.  At  time  t = 0,  the  illumination  ends,  and  the 
generation  rate  vanishes.  Then  (4.1)  becomes 


5n 

at 

n 

Tn 

+ D v2n 
n 

(4.2) 

For  the 

sake 
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sin  exp 


(-Mijkt)  (4.3) 


where 


(4.4) 
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anc*  Gi  jk  amPl itude  of  the  ijk  mode.  The  reciprocal  of  M-  .. 

defines  the  lifetime  of  the  excess  carrier  density  for  a given  mode 
(ijk)  in  the  filament.  In  an  investigation  similar  to  that  of  [42], 
Shockley  [6]  argued  that  the  lifetime  of  the  lowest  mode  is  much  longer 
than  those  of  the  others;  all  terms  in  (4.3)  will  die  away  rapidly 
except  the  first  one.  Stevenson  and  Keyes  [42]  found  that  the  lifetime 
of  some  modes  are  comparable  to  that  of  the  first  mode,  but  that  the 
amplitudes  of  the  higher  modes  are  much  smaller  than  that  of  the  first 
one.  Thus,  the  lifetime  of  the  lowest  mode  is  usually  defined  as  the 
lifetime  of  the  filament  under  the  boundary  conditions  defined 
previously: 


it2D 


xf  ~ Mlll 


n /I 


• tr + V1  (-T  * T + 

Tn  4 A2  B2  C2 


(4.5) 


4.2.2  Effective  Filament  Lifetime  at  Steady  State 

We  now  try  to  use  the  first-term  approximation  to  define  the 
effective  filament  lifetime  at  steady  state.  Let  the  filament  be  the 
quasi -neutral  base  of  a forward  biased  p/n  junction;  the  face  z = 0 is 
the  base  edge  of  the  p/n  junction  space-charge  region.  Then  (4.2) 
becomes 


0 = --  +Dv2n 
n 


(4.6) 


with  boundary  conditions: 
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(4.7) 


at  the  injecting  face  z - 0 and  S -*•  <*>  on  the  other  five  faces  of  the 
filament.  The  solution  is 


n(x,y,z)  = 


n4 


i odd  j odd  NAA(0)  ^ TJ  sin  sin  ^ 


sinh[ri j (2C-z) ] 


2B  sinh  2Cr. . 

TJ 


(4.8) 
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(4.9) 


By  the  first-term  approximation,  (4.8)  becomes 


n(x,y,z)  = nu(x,y,z) 


n2 


xy  s1nh[ru(2C-z)] 
2B  sinh  20^ 


(4.10) 


The  current  density  on  the  injecting  face  is  given  by 
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J(x,y,0) 


= qD 


an 
n 5z 


(4.11) 


and  the  total  current  in  the  base  region  has  the  form 
2A  2B 

K=0  = / / J(x,y,0)  dxdy  . (4.12) 

z 0 0 

Combining  (4.10),  (4.11),  (4.12)  gives 

1/2 

(4.13) 
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The  mathematical  form  for  the  total  current  is 


z=0 


,1/2 

n 

1/2 


(4.14) 


as  suggested  in  [45],  Hence  the  effective  filament  lifetime  in  the 
first-term  approximation  can  be  defined  as 


= [M 


11 


2CM^ 

coth2  i-jif)]-1 


(4.15) 


Note  that  (4.5)  and  (4.15)  are  equal  for  a semi -i nfi nite  filament,  and 
that  Stevenson  and  Keyes  [42]  have  suggested  the  adequacy  of  (4.5)  for 
characterizing  filiment  transient  response.  In  the  semi-infinite  case 


(C  -►  «),  the  current  density  Jz=0  in  a one-dimensional  diode  (A  -►  =°, 
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B »)  for  the  steady-state,  obtained  from  (4.13)  or  from  (4.14)  and 
(4.15),  is 


J 


z=0 


lz=0  _ 64  r qni 
4AB  ' ^ { N^TOT 


[exp(^y)  - 1] 


)1/2 
1/  t) 


(4.16) 


where  the  factor  enclosed  by  curly  brackets  is  recognized  as  the 
customary  one-dimensional  current  density.  The  disagreement  of  (4.16) 
by  a multiplicative  factor,  64/ it4,  for  the  steady-state,  prompts  our 
further  assessment  of  the  quality  of  the  first-term  approximation,  which 
now  follows. 


4.2.3  Divergence  of  Exact  Solution 

To  examine  the  quality  of  the  first-term  approximation,  consider 
the  exact  solution  yielded  by  combining  (4.8),  (4.11),  and  (4.12).  A 
method  to  sum  the  series  is  not  apparent.  Thus  we  deal  with  the  issue 
of  term-by-term  differentiation  of  the  series,  integrating  the  result  to 
determine  I(V). 

To  simplify  the  discussion,  without  loss  in  generality,  we  treat  a 
simpler,  related  boundary-value  problem  yielded  by  letting  % become 
infinite  and  letting  the  structure  become  semi-infi nite  in  the  x and  z 
directions  while  requiring  that  n(x,y,0)  on  the  injecting  face  be 
independent  of  x and  y and  n = 0 on  the  two  lateral  sides.  The  result 
is  a boundary-value  problem  for  which  n(y,z)  satisfies  Laplace's 
equation  in  two  dimensions.  Analogous  problems  exist:  for  example,  a 

semi -i nfi nite  rectangular  prism  of  thermally  (electrically)  conducting 
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material  in  which  the  top  side  is  held  at  constant  temperature 
(potential)  and  the  lateral  sides  are  maintained  at  zero  temperature 
(potential)  [41,46],  The  solution  is 


"(y’z)  = j odd 


n? 


j odd  Naa(0)  texP(kT^  ^ u j exp  S1n  (^B^ 


(4.17) 


which  can  be  summed  to  yield  a closed-form  solution  [41,46],  also 
obtainable  by  comformal  mapping  [47]: 
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(4.18) 


Differentiating  (4.18)  with  respect  to  z at  z = 0,  we  obtain 
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(4.19) 


There  are  two  singularities  near  (y=0,  z=0)  and  (y=2B,  z=0)  in  (4.19). 
In  the  semiconductor  problem,  integrating  (4.19)  along  the  base  edge  of 
the  p/n  junction  space  charge  region  from  y = 0 to  y = 2B,  we  obtain  an 
infinite  total  current. 

Based  on  a series  like  (4.17),  one  might  think  that  the  total 
current  converges.  Three  erroneous  lines  of  reasoning  might  lead  to 
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this  conclusion.  First,  one  could  assume  that  (4.17),  or  its  analog  in 
the  time-invariant  semiconductor  filament  problem,  can  be  differentiated 
termwise  with  respect  to  z.  But  inspection  of  (4.17)  shows  that  each 
term  of  the  differentiated  series  vanishes  at  (y=0,  z=0)  and  (y=2B, 
z=0),  in  conflict  with  the  singularities  at  these  two  points  implied  by 
(4.19).  This  invalidates  termwise  differentiation.  Second,  one  could 
investigate  numerically  the  sequence  of  partial  sums  differentiated  with 
respect  to  z.  This  differentiated  sequence  may  appear  to  converge  at 
(0,0)  and  (2B,0)  if  one  takes  too  few  terms,  resulting  in  the  incorrect 
conclusion  that  the  current  I(V)  converges.  Third,  as  in  [7],  one  could 
seek  a numerical  solution  by  finite-difference  techniques  of  the  model 
under  study  in  which  infinite  S exists  on  the  lateral  and  back 
surfaces.  The  exact  solution  of  this  model  diverges,  as  we  have 
shown.  Thus,  only  by  some  numerical  subtlety  could  one  obtain  a 
convergent  solution,  and  the  physical  meaning  of  such  a solution,  as 
obtained  in  [7],  for  3,000  grid  points,  is  unclear. 

The  merit  of  examining  (4 .17) -( 4.19)  in  the  simpler,  related 
problem  comes  from  the  ease  in  identifying  the  singularities  at  (0,0) 
and  (2B,0)  as  the  origin  of  the  divergences  of  the  I(V)  characteristic. 
These  singularities  remain  in  the  three-dimensional  semiconductor 
filament  problem  because  introduction  of  finite  lifetime  and  increase  of 
space  dimensionality  does  not  alter  the  discontinuity  in  n at  the  edges 
of  the  injecting  face  in  the  model  used  (S  ® on  lateral  surfaces). 
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4.2.4  Fourier-Series  Approximation 

The  divergence  just  discussed  comes  from  a flaw  in  the  model 
used.  Often  in  semiconductor  device  problems,  one  assumes  infinite 
surface  recombination  velocity,  implying  zero  excess  carrier  density,  to 
get  a simple  result.  An  example  is  the  short-base  p/n  junction  diode  in 
which  one  assumes  that  the  ohmic  contact  at  the  back  surface  maintains 
the  carrier  densities  there  at  thermal-equilibrium  values. 

But  whatever  mechanisms  at  an  ohmic  contact  act  to  reduce  the 
excess  carrier  densities,  the  result  is  never  reduced  to  zero.  For 
example,  the  unilateral  thermal  velocity  is  finite,  recombination  rates 
are  finite,  and  quantum-mechanical  tunneling  never  gives  absolute 
removal . 

For  similar  reasons,  the  surface  recombination  velocity  on  the 
lateral  surfaces  of  the  semiconductor  filament  always  must  remain 
finite.  Thus  the  discontinuity  in  n at  the  edges  of  the  injecting  face, 
present  up  to  now  in  the  model  we  and  others  have  used,  does  not 
physically  exist,  and  the  resulting  singularities  in  dn/dz  and  the 
resulting  divergence  of  the  I(V)  characteristic  comes  from  the 
inappropriateness  of  the  model  used. 

The  most  straightforward  way  of  modifying  the  model  in  the  time- 
invariant  semiconductor  filament  problem  consists  in  making  S finite. 
But  this  yields  a boundary-value  problem  for  which  a closed-form 
solution  has  not  appeared,  requiring  therefore  treatment  by  numerical 


methods . 
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We  now  develop  a modification  of  the  model  that  leads  to  simple, 
approximate  solutions  in  analytic  form.  The  excess  electron  density  at 
z=0,  0 < y < 2B,  can  be  represented  in  terms  of  Fourier  sine  series  [48] 
as  follows: 


Let  any  partial  sum  Sj  of  (4.20)  be  taken  as  an  approximation;  Sj 

vanishes  at  z=0,  y=0  and  z=0,  y=2B,  and  no  discontinuities  result  if  we 

continue  to  use  the  approximation  S -*•  °°  on  the  lateral  surfaces.  In 

this  model,  for  xn  -►  ® and  a semi -infinite  region  in  the  z direction, 

the  solution  is  (4.17)  at  term  J because  of  the  orthogonality  property 
of  the  sine  function.  The  truncated  solution  satisfies  the  modified 
boundary  condition  in  the  closed  interval  z=0,  0 < y < 2B,  and  can  be 
differentiated  or  integrated.  At  the  mid-point,  y=B,  if  the  average  of 
any  two  consecutive  partial  sum  is  taken,  the  boundary  value  of  n on  the 
injecting  face  is  nearly  the  left  side  of  (4.20)  [49].  Thus  a simple 
form  of  the  truncated  Fourier  sine  series  of  the  boundary  condition  at 
z=0  is  the  average  of  the  second  and  third  partial  sums  as  follows: 


(4.20) 


TI^Tor  [exP(kT>  - 1](;  <s1n 


(4.21) 


This  truncated  Fourier  sine  series  is  shown  in  Fig.  4.1. 
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4-2.5  Definition  of  Filament  Lifetime  for  Steady  State  and  S Infinite 
on  Lateral  Surfaces 

To  simplify  the  discussion  we  continue  treating  a two-dimensional 
model.  Extension  to  three  space  dimensions  is  straightforward. 

We  now  include  finite  xn,  so  that  the  Laplace's  equation  treated 
earlier  in  Sec.  4.2.3  becomes  generalized  to 


a2n  52n  n 
by2  5z2  ^nTn 


(4.22) 


The  key  to  the  method  is  the  use  of  the  Fourier-series  approximation  of 
Sec.  4.2.4  to  give  for  the  boundary  conditions  on  the  injecting  face. 


n = 


AA 


Toy  ^exP(kT") 


1]  1 [sin  2B  ' 3 2B 


4 sin  4rcy  + 

3 2B  2x5 


sin  2^-j  at  z~ 


n = 0 at  y = 0 and  y = 2B 


(4.23) 

(4.24) 


For  the  average  of  partial  sums  illustrated  in  Figure  4.1,  the  solution 
for  (4.22)  and  (4.23)  is 
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Fig.  4.1  The  first  partial  sum  and  two  averages  of  partial  sums 
compared  with  the  space-independent  (normalized)  boundary 
condition  on  the  injecting  face.  The  highest-order  average 
of  partial  sums  shown  approximates  the  space-independent 
boundary  condition  well  (within  3%)  over  about  half  of  the 
face  near  the  mid-point,  and  merges  wihtout  discontinuity 
with  the  boundary  conditions  on  the  lateral  surfaces,  but  it 
exhibits  Gibbs  phenomenon  resulting  in  an  error  of  15%  at 
the  two  maxima. 
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where  M . = f—  + 

1 ' *r 


(4.26) 


Thus  the  total  current  is 


(4.27) 


where 


(4.28) 


defines  the  filament  lifetime  for  steady-state  conditions,  for  a model 
in  which  S + - on  the  lateral  surfaces.  Here  we  use  the  same  symbol  as 
that  in  (4.5).  Equation  (4.28)  is  an  example  illustrating  the  new 
definition  of  the  steady-state  filament  lifetime  for  a partial-sum 
approximation  of  the  Fourier  series  of  (4.20);  it  involves  partial  sums 
including  the  third  term  of  the  Fourier  series.  To  develop  expressions 
for  the  filament  lifetime  corresponding  to  additional  higher-order  terms 
of  this  series,  one  can  follow  the  procedure  outlined  in  Sec.  4.2.4;  the 
inclusion  of  additional  higher-order  terms  does  not  necessarily  imply 

better  accuracy  because  of  the  Gibbs  phenomenon  [50]  displayed  in  Fig. 

4.1. 

It  is  instructive  to  compare  the  steady-state  filament  lifetime  of 
(4.28)  with  the  transient  filament  lifetime  of  Shockley  [6],  which 
Fig.  4.2  illustrates.  Note  the  agreement.  A main  purpose  of  the 
foregoing  development  has  been  to  demonstrate  that  Shockley's  filament 
lifetime  holds  as  an  approximation  for  the  time-invariant  steady  state. 


Fig.  4.2  Comparison  of  the  Shockley  filament  lifetime  [6]  with 
the  filament  lifetime  of  Eq.  (4.28)  of  the  present 
work.  Low  injection  in  a p-type  semiconductor  is 
assumed  for  a filament  having  infinite  surface 
recombination  velocity  on  the  lateral  surfaces.  The 
lifetime  of  excess  electrons  appears  as  a parameter; 
the  diffusivity  Dn  = 25  cirr/s  in  the  calculations. 


(O0S)  *_L 
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a conclusion  not  necessarily  to  be  expected  from  the  reasoning  used  by 
Shockley  [6]  and  others  [42-44].  This  is  an  important  conclusion, 
partly  because  other  workers  [45,  51]  have  assumed  this  is  so;  the 
present  work  provides  a basis  for  this  assumption. 

4.3.  Discussion 

For  a filament  having  S=»  on  the  lateral  surfaces,  the  first-term 
approximation  suggested  by  Shockley  is  meant  to  define  the  effective 
lifetime  only  for  transient  conditions.  For  the  time-invariant  case 
(steady  state),  we  have  shown  that  the  infinite-series  form  of  solution 
for  the  excess  minority-carrier  distribution  (in  a forward-biased  p/n 
junction  solar  cell,  for  example)  results  in  infinite  current.  The 
demonstration  here  treats  the  steady-state  injection  of  excess  carriers 
at  a front  surface  such  as  that  adjoining  a forward-biased  p/n  or 
Schottky-barrier  junction.  The  injected  carriers  recombine  in  the 
volume  or  at  the  lateral  surfaces.  Such  problems  arise,  for  example,  in 
the  analysis  of  polycrystall i ne  solar  cells,  for  which  the  lateral 
surfaces  account  for  the  presence  of  grain  boundaries.  Consideration  of 
such  problems  has  given  rise  to  the  method  proposed  here  involving  a 
Fouri er-seri es  approximation  of  the  space-independent  boundary  condition 
on  the  injecting  face. 

The  approach  used  in  this  study  resembles  an  approach  used  by 
Kantorovich  and  Krylov  [52].  These  workers  discuss  various  methods  for 
speeding  the  convergence  of  series  that  result  from  solving  boundary- 
value  problems  related  to  the  one  we  study  here.  The  general  approach 
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used  involves  a prel imi nary  improvement  of  the  boundary-value  problem  at 
the  outset  that  will  assure  rapid  convergence.  The  aim  of  the 
improvement  is  to  change  the  boundary-value  problem  to  yield  existence 
and  continuity  of  the  first  derivatives  of  (a)  the  boundary  function  and 
(b)  the  inhomogeneous  part  of  the  differential  equation.  They  give  a 
general  method  for  accomplishing  this.  Our  approach  is  the  same  as 
theirs  in  the  sense  that  we  change  the  boundary  function  to  yield 

existence  and  continuity  of  the  first  derivatives.  The  change  we  make 

removes  delta  functions  in  the  first  derivative  on  the  boundary.  We 

emphasize  that  this  change  alters  the  mathematical  model  of  the  physical 
system  we  treat  in  a direction  toward  the  actual  physical  system.  The 
actual  physical  system  does  not  contain  singularities  at  the  end  points 
of  the  boundary  because  infinite  surface  recombination  velocity  is  not 
physically  realizable. 

Nonetheless,  the  treatment  here  of  a model  having  infinite  surface 
recombination  velocity  on  the  lateral  surfaces  yields  several  benefits. 
First,  we  examined  a model  used  by  previous  workers,  who  did  not 

apparently  inquire  into  the  convergence  of  the  exact  solution.  Second, 
for  the  steady-state  condition  treated,  we  demonstrated  that  for  such  a 
model  the  exact  solution  diverges,  indicating  the  advisability  of 
caution.  Third,  through  the  heuristic  approximation  made  in  treating 
this  model,  we  have  enabled  the  condition,  S = »,  to  be  preserved  while 
obtaining  a closed-form  convergent  solution;  among  other  benefits,  this 
has  the  virtue  of  introducing  zero-density  conditions  on  the  lateral 
surfaces,  which  lend  themselves  to  the  use  of  the  superposition 
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principle  in  constructing  solutions  to  boundary-value  problems  of 
interest  in  semiconductor  device  design,  as  in  Sec.  4.4.3  that  follows. 

4.4  Ap  pi i cations 

This  section  indicates  various  uses  of  the  method  described  here 
for  a variety  of  problems.  In  all  of  these  problems,  the  filament 
lifetime  follows  straightforwardly  from  applying  the  procedure  of 
Sec.  4.2.5. 


4-4*l  First  Example:  Three-Dimensional  Rectangular,  Finite  Filament 

The  method  of  Sec.  4.2  is  easily  extended  to  treat  a three- 
dimensional  boundary-value  problem.  Consider  the  filament  illustrated 
in  Fig.  4.3.  The  injecting  face  is  the  plane  defined  below  in  (4.29), 
for  which  the  method  of  Sec.  4.2  yields  a Fourier-series  approximation 
of  the  boundary  condition  given  by 
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0 < y < 2B,  0 < x < 2A,  z=0 


(4.29) 
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Fig.  4.3  The  boundary  conditions  on  the  surfaces 
of  a rectangular  parallel  piped  filament. 
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The  filament  lifetime  then  follows  straightforwardly  from  the  procedure 
of  Sec.  4.2.5.  This  is  the  filament  lifetime  of  a nonil  luminated 
rectangular  grain  for  which  the  potential  barrier  at  the  face  indicated 
in  (4.29)  is  lowered  by  the  applied  forward  bias  V and  for  which  the 
surface  recombination  velocity  on  the  lateral  surfaces  is  approximately 
infinite.  The  boundary  condition  on  the  plane  corresponding  to  z = 2C 
is  appropriate  for  an  ohmic  contact  at  that  plane. 


Second  Example:  Cylindrical,  Finite  Filament  with  Ohmic  Contact 
on  the  Back  Surface 

The  approximation  for  the  boundary  condition  on  the  injecting 
surface  shown  in  Fig.  4.4  is  a Bessel -function  polynomial: 


r0r01Jl(r01r0)  r0r02Jl(r02r0) 


+ Vr03rQ)  ^ 


(4.30) 


where  JQ  and  are  the  zeroth-order  and  the  first-order  Bessel  function 
respectively,  and  rQm  are  the  r that  make  JQ(rr)  = 0. 

The  solution  for  original  boundary  condition  (Fig.  4.4)  is  the 
infinite  series 
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Fig.  4.4 


The  boundary  conditions  on  the  surfaces 
of  a cylindrical  filament. 
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The  partial  sum  in  (4.30)  is  justified  for  the  same  reason  as  that 
indicated  between  (4.20)  and  (4.21)  in  Sec.  4.2.4.  Figure  4.5  shows  the 
original  and  the  modified  boundary  functions.  This  example  is  the 

counterpart,  for  a cylindrical  grain,  of  the  first  example  in 
Sec.  4.4.1,  subject  to  the  assumption  that  n has  no  angular  dependence. 
The  model  treated  here  is  the  same  as  that  treated  in  [7]  by  standard 
numerical  finite-difference  techniques. 

4*4*3  Third  Example:  Superposition  of  Elemental  Solutions  to  Solve 
a Class  of  Device  Problems 

Consider,  for  example,  a cylindrical  grain  of  p-type  semiconductor. 
Suppose  that  the  phosphorus  diffusion  that  forms  the  n/p  junction,  and 
the  injecting  face  for  the  quasi-neutral  grain,  also  provokes  a 
preferential  diffusion  along  the  grain-boundary  surfaces  [53].  This 
preferential  diffusion  forms  an  n/p  injecting  face  along  the  lateral 
surfaces  nearly  extending  to  the  back  surface.  Suppose  further  that  the 
region  near  the  back  surface  has  an  effective  recombination  velocity  of 
nearly  zero.  This  can  be  achieved  by  a low-high  junction  (p/p+)  as  in  a 
back -surface-field  solar  cell  or  by  other  means.  The  resulting  device 
can  be  modeled  as  shown  in  Fig.  4.6(a).  As  a solar  cell,  the  non- 
ill  uminated  current  can  be  small  if  the  electron  diffusion  length  Ln 
(assuming  low  injection  conditions)  is  much  less  than  both  rQ  and  2C. 

Exploiting  the  linear  relation  between  n and  Jn  to  employ  the 
superposition  principle,  we  determine  the  solution  for  the  problem  of 
Fig.  4.6(a)  by  adding  the  solutions  for  the  full -line  portions. 
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Uf) 


The  first  partial  sum  and  two  averages  of  partial  sums  compared  with  the  space- 
independent  (normalized)  boundary  condition  on  the  injecting  face  for  a cylindrical 
grain. 
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excluding  the  dashed-line  portion  of  Figs.  4.6(b)  and  4.6(c).  The 
condition  S = 0 at  z = 2C  results  from  symmetry  in  the  solution  of 
Figs.  4.6(b)  and  4.6(c)  if  the  dashed-line  portions  are  included.  For 
Fig.  4.6(b)  we  use  a truncated  Bessel -series  approximation,  as  in 
Sec.  4.4.2  and  Fig.  4.5,  and  for  Fig.  4.6(c)  we  use  a truncated  Fourier- 
series  approximation. 

If  in  the  device  problem  of  Fig.  4.6(a),  we  were  to  replace  the 
ideal  low-high  junction  by  an  ohmic  contact,  then  the  solution  results 
from  superposing  solutions  for  Fig.  4.4  and  Fig.  4.6(c)  (including  the 
dashed-line  portion,  with  appropriate  scaling). 

Rectangular  geometries  and  other  related  problems  yield  to  similar 
superpositions  of  elemental  solutions. 


Fig.  4.6 


The  boundary  conditions  on  (a)  are  equivalent  to  the 
superposition  of  the  full -line  portions  of  (b)  and  (c). 
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CHAPTER  5 

FINITE  SURFACE  RECOMBINATION  VELOCITY  AT  GRAIN  BOUNDARIES 

5.1  Introduction 

We  pointed  out  in  Chapter  4 that  the  assumption  of  infinite  surface 
recombination  velocities  on  the  lateral  surfaces  of  the  semiconductor 
filament  results  in  divergence  of  the  I ( V)  characteristic.  Furthermore 
we  showed  that  the  introduction  of  truncated  Fourier  sine  or  Bessel  - 
function  series  to  modify  the  original  boundary  conditions  eliminates 
the  singularities  in  differentiation  (dn/dz)  and  gives  analytic 
solutions  for  mi nority-carrier  density,  total  current,  and  effective 
filament  lifetime  in  a semiconductor  filament. 

The  Fourier-series  (Bessel -function)  approximation  suggested  in 
Chapter  4 applies  strictly  for  S = ® and  has  been  assumed  adequate  for  a 
semiconductor  filament,  with  a typical  high  value  of  surface 
recombination  velocity  ( 106  cm/s)  on  the  lateral  surfaces.  Any  process 
making  the  surface  recombination  velocity  much  smaller  than  106  cm/s 
invalidates  the  Fourier-series  (Bessel-function)  approximation 
of  Chapter  4. 

Thus  methods  of  solving  the  boundary-value  problem  for  arbitrary 
surface  recombination  velocity  on  the  lateral  surface  are  of  interest. 
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They  give  the  I(V)  characteristic  for  small  surface  recombination 
velocity  (<  106  cm/s),  which  cannot  be  treated  by  the  method  of 

Chapter  4,  and  they  offer  a comparison  between  the  results  of  Fourier- 
series  (Bessel -function)  approximation  and  those  when  the  finite  surface 
recombination  velocity  is  set  to  be  106  cm/s. 

Efforts  have  been  made  to  solve  the  filament  boundry-val ue  problem 
with  arbitrary  but  position  independent  surface  recombination  velocity 
on  the  lateral  surfaces. 

We  will  focus  on  a method  due  to  Gonzalez  [8].  Gonzalez 
demonstrated  that  the  boundary-value  problem  mentioned  above  is  just  the 
Sturm-Liouvill e problem  [54].  Following  standard  procedures  for  the 
Sturm-Liouville  problem,  he  obtained  an  infinite-series  solution  for  the 
mi nority-carrier  concentration.  Expecting  finite  total  current  at  the 
injecting  surface  under  the  assumption  of  finite  surface  recombination 
velocity  on  the  lateral  surface,  he  first  differentiated  term-by-term 
the  infinite  series  describing  the  minority  carrier  concentration  to 
obtain  current  density.  Then  he  integrated  term-by-term  the  infinite 
series  describing  the  current  density  to  obtain  the  total  current. 
Gonzalez  did  this,  however,  without  checking  the  uniform  convergence  of 
the  infinite  series,  which  is  a necessary  procedure  [55]  but  apparently 
a difficult  one. 

The  purpose  of  this  chapter  is  to  indicate  the  validity  of  term-by- 
term  differentiation  and  integration  while  avoiding  the  difficult  issue 
of  uniform  convergence.  We  base  our  reasoning  on  an  implication  of  the 
development  of  Fourier-series  approximation  of  Chapter  4. 
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5.2  Sturm-Li ou vi 1 1 e Problem 

Consider  the  semi -infinite  quasi-neutral  base  of  a forward-biased 
two-dimensional  p/n  junction  diode;  the  differential  equation  governing 
the  excess  minority  electron  is 


= 0 


(5.1) 


with  boundary  conditions 


n(y,0)  = Cexp(kT}  ’ 1] 


on  the  injecting  size  z = 0 and 
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on  the  two  lateral  sides,  as  shown  in  Fig.  5.1. 

Following  the  procedure  for  solving  the  Sturm-Liouville  problem 
[54,56],  we  solve  (5.1)  with  boundary  conditions  (5.2),  (5.3),  and  (5.4) 
and  obtain  an  infinite-series  solution  for  minority-carrier  density  in 
quasi-neutral  base  of  a forward-biased  two-dimensional  p/n  junction 
diode  as  follows: 
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Fig.  5.1  The  boundary  conditions  for  the  semi -infinite  quasi-neutral 
base  of  a forward-biased  two-dimensional  p/n  junction  diode. 
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and  F-  can  be  obtained  by  solving 
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A method  to  sum  a complicated  series  like  (5.5)  is  apparently  not 
available.  This  difficulty  led  Gonzalez  to  differentiate  the  series 
term  by  term  and  then  integrate  the  result.  This  gives  rise  to  the 
issues  of  mathematical  rigor  discussed  above. 

We  now  review  the  boundary-value  problem  described  in  Sec.  4.2.3. 

The  relevance  of  this  review  to  the  problem  now  considered  will  become 
apparent. 

This  problem  (Sec.  4.2.3)  is  basically  the  same  as  the  problem  we 
now  treat  except  that  the  recombination  inside  the  boundary  was 
neglected  in  Chapter  4 and  the  surface  recombination  velocity  on  the 
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lateral  sides  was  assumed  to  be  infinite.  In  fact,  the  method  employed 
to  yield  the  solution  (4.17)  for  infinite  S is  also  the  same  procedure 
used  in  Sturm-Liousvil  le  (S-L)  problem,  because  infinite  surface 
recombination  implying  zero  concentration  of  excess  minority  on  the 
lateral  sides  is  a special  case  for  boundary  conditions  in  the  S-L 

problem  [54,  56],  The  purpose  of  pointing  out  the  utilization  of  the 

same  method  in  these  two  problems  is  to  try  to  find  some  way  of 
assessing  the  convergence  of  the  infinite-series  solution  for  finite  S 
boundary  conditions  from  the  experience  gained  in  determining  the 
divergence  of  exact  solution  for  infinite  S on  the  lateral  sides 

(Chapter  4) . 

We  have  demonstrated  that  the  solution  for  infinite  S is  also  an 
infinite  series,  and  that  it  can  be  summed  to  yield  a closed-form 

solution  in  which  the  existence  of  discontinuities  near  the  corners  of 
the  injecting  side  results  in  infinite  current  density.  As  we  mentioned 
earlier,  a closed-form  solution  is  not  present  for  finite  S.  Hence  we 

have  no  opportunity  to  assess  the  convergence  of  the  solution  by  this 
way. 

The  divergence  of  the  solution  for  infinite  S can  also  be 
determined  by  examining  the  infinite-series  solution.  Letting  z = 0, 
y - 0 (or  2B)  in  (4.17),  we  find  that  each  term  of  the  series  is  zero. 
Hence  the  excess 2minority  carrier  density  is  zero  at  (y=0,  z = 0),  but 
is  constant  ^N'A^(o)^exP^kT^  " ^ ln  the  interval  0 < y < 2B  on  the 
injecting  side,  because  the  method  employed  to  solve  the  problem  ensures 
the  solution  satisfying  the  boundary  conditions  in  that  interval. 
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The  abrupt  change  on  excess  minority  carrier  density  near  the  corners  on 
the  injecting  side  is  the  origin  of  the  discontinuities. 

Using  the  approach  just  described,  we  let  z = 0,  y = 0 (or  2B)  in 
infinite  series  (5.5)  and  have  a series  with  simpler  form  in  every 
term.  Again  the  new  series  is  not  easily  summed.  Not  knowing  the 
excess  minority-density  at  the  corners  of  the  injecting  side,  we  are  not 
able  to  make  any  conclusion  about  the  continuity  of  the  solution  in  the 
injecting  side. 

From  Chapter  4,  however,  it  is  obvious  that  the  discontinuities  for 
infinite  S occur  in  one  spatial  dimension,  on  the  injecting  side  z = 0, 
although  the  original  problem  is  in  two-dimensional  space.  This  gives 
us  a hint  that  we  might  examine  the  divergence  of  the  corners  by  using  a 
one-dimensional  model  as  follows. 

Consider  the  quasi-neutral  base  of  a forward  biased  one-dimensional 
p/n  junction  diode  whose  back  is  covered  with  an  ohmic  contact. 
Assuming  surface  recombination  velocity  on  the  back  surface  is  infinite, 
we  have 

n2 

nU)  = VaA(0)  CeXP(kT)  " 1](WB  ' x)  (5.9) 


for  the  excess  minority-carrier  density  in  a short  base  of  thickness 
Wg.  The  electron  density  at  the  base  edge  of  the  space-charge  region 
x = 0 is 

n(0)  = N^Wexp(kT  ' 1)]  ’ 


(5.10) 
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and  that  on  the  back  surface  x = WD  is 

D 


n(WB)  = 0 


(5.11) 


If  the  base  width  shrinks  to  zero,  a case  that  the  ohmic  contact  is 
directly  in  touch  with  the  base  edge  of  the  space-charge  region,  (5.11) 
becomes 


n(0)  = 0 


(5.12) 


The  contradiction  between  (5.10)  and  (5.12)  implies  that  electron 
density  is  a step  junction  at  x = 0. 

In  the  original  two-dimensional  problem,  the  base  edge  of  the 
space-charge  region  contacts  the  surface  with  infinite  S at  the  corners 
(y  = 0,  z = 0)  and  (y  = 2B,  z = 0) . Using  the  model  just  developed,  we 
found  that  discontinuities  arise  at  those  two  points,  which  agrees  with 
the  conclusions  obtained  by  previous  two.  methods. 

Now  consider  the  same  one-dimensional  problem  again  with  finite  S 

on  the  back  surface.  The  minority  electron  density  in  the  short  base  is 
then 
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On  the  back  surface. 
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(5.14) 


Letting  Wg  be  zero,  we  have  a constant  (yj— joy[exp(^)  - 1)])  electron 
density,  the  same  result  as  that  obtained  by  substituting  x = 0 into 
(5.13).  This  is  a quite  different  result  from  that  under  the  assumption 
of  infinite  S.  The  simplicity  of  the  one-dimensional  model  makes 
evident  that  S < °°  provides  a transition  to  a convergent  solution. 

Using  this  model  as  a guide,  consider  the  electron-density  function 
at  the  corners  of  the  injecting  side  in  the  two-dimensional  problem  in 
which  S < °°.  The  electron  density  is  constant  for  0 < y < 2B.  Because 
S < ”,  by  analogy  we  see  that  no  discontinuities  exist  at  the  corners  of 
the  injecting  side.  Thus  the  analogy  strongly  suggests  that  termwise 
differentiation  of  the  infinite  series  (5.5)  and  termwise  integration  of 
the  result  will  not  cause  any  mathematical  difficulty.  Thus  the  finite 
total  current  is 
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where  Gj , Aj , and  V-  has  been  defined  by  (5.6),  (5.7)  and  (5.8) 

respectively.  The  effective  filament  lifetime  can  be  defined  by  charge- 
control  model  as 
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where 
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5 .3  Comparisons 

Table  5.1  lists  I and  xf  calculated  by  finite-S  method  described  in 
the  foregoing  sections  and  by  Fourier-series  approximation,  in  which  S 
is  assumed  to  be  infinite  (Chapter  4)  for  different  values  of  2B.  Where 
Naa  and  finite  S are  assumed  to  be  1017  cm-3  and  106  cm/sec 
respectively,  Dp  is  calculated  by 


(kT/q){92  + (1360  - 92)/[l  + (NAA/1.3  x 1017)0-91]}  [21], 

and  xn  has  the  value  of  10“3  sec  [18]. 

We  have  indicated  the  convergence  of  I for  assumption  of  finte  S. 
The  accuracy  of  the  result  increases  when  more  terms  are  summed.  Here 
we  take  the  600th  partial  sum  for  comparison. 

In  Table  5.1,  we  see  that  Fourier-series  is  a good  approximation 
for  both  effective  filament  lifetime  and  total  current,  especially  when 
a quick  estimation  is  required,  because  of  its  simplicity. 


5 .4  Discussion 

The  main  contribution  of  this  chapter  is  the  demonstration  that  the 
method  of  Chapter  4 for  S = » provides  a good  approximation  for 
filaments  for  which  S ~ 106  cm/s.  This  demonstration  followed  from 
consideration  of  an  analogous  one-dimensional  problem  in  which  S < <=  or 
S = 00  played  a prominent  role. 
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Table  5.3 

Illustrating  the  adequacy  of  the  approximation  of  Chapter  4, 
which  improves  as  the  width  2B  increases  and  reduces 
the  influence  of  the  corners 
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CHAPTER  6 

SUMMARY  AND  DISCUSSION 

There  is  interest  in  designing  Si  solar  cells  made  from  inexpensive 
processing  and  materials.  Such  solar  cells  could  have  a high 
concentration  of  recombination  centers  throughout  the  cell.  In  this 
dissertation  we  explored,  through  both  experimental  and  theoretical 
study,  several  design  principles  that  have  been  previously  suggested 
[35,57]  for  obtaining  high  power  conversion  efficiency  from  p/n  junction 
Si  solar  cells  that  contain  a high  concentration  of  recombination 
centers  throughout  the  active  base  region.  These  design  principles 
include  (a)  an  n-type  base  region  having  a graded  dopant-atom  profile 
that  increases  from  the  front  surface  toward  the  back  surface,  yielding 
a drift-field  (DF)  solar  cell;  (b)  cell  thinness;  (c)  gettering  of 
metallic  and  other  impurities  to  the  n+  region  bordering  the  back 
surface  of  the  cell;  (d)  treatment  of  grain  boundaries  by  preventing  the 
formation  of  p/n  junctions  at  grain-boundary  surfaces,  or  by  exploiting 
the  rapid  diffusion  of  donor  or  acceptor  atoms  down  grain  boundaries. 

We  separated  the  description  of  the  research  into  two  parts. 
First,  we  evaluated  the  performance  of  single-crystal  Si  DF  solar 
cells.  Chapter  2 dealt  with  the  theoretical  portion  and  Chapter  3 with 
the  experimental  portion  of  the  work.  Second,  we  investigated  the  I-V 
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characteristics  of  polycrystall  i ne  Si  solar  cells  (without  DF).  In 
Chapter  4 we  obtained  closed-form  solutions  for  the  I-V  characteristics 
by  combining  the  assumption  of  infinite  surface  recombination  velocity  S 
on  the  grain  boundaries  with  a modified  boundary  condition  at  space- 
charge-region  edge  of  the  p/n  junction.  In  Chapter  5,  for  comparison, 
we  treated  a device  having  finite  S.  Here  S is  the  effective  surface 
recombination  velocity  at  the  space-charge-region  edge  of  the  grain 
boundary. 

In  Chapter  2,  we  presented  general  analytic  solutions  for  static 
I-V  characteristics  of  quasi-neutral  regions  of  non-il luminated  Si  DF 
solar  cells  under  the  following  assumptions:  (a)  the  quasi-neutral 
region  has  a graded  shallow-level  impurity  concentration  producing  a 
constant  built-in  electric  (drift)  field;  (b)  minority  carriers  injected 
into  this  region  stay  at  concentrations  low  enough  to  avoid  violation  of 
low-injection  conditions;  (c)  the  mi  nority-carrier  lifetime  of  this 
region  depends  on  position  in  accordance  with  a power-law  dependence  on 
the  shallow-level  donor  concentration,  a dependence  that  is  consistent 
with  the  longest  minority-carrier  lifetimes  measured  and  with  the 
physical  chemistry  of  di vacancy-donor  reactions  at  high  temperatures. 
The  solutions  presented  are  apparently  the  first  that  include  assumption 
(c)  above.  The  solution  showed  that  (a)  the  pre-exponential  factor  jpo 
for  the  hole  current  density  in  the  non-il luminated  cell  under  forward 
voltage  V is  insensitive  to  the  surface  recombination  velocity  S on  the 
back  surface  (this  implies  that  no  additional  processes  are  needed  for 
reducing  S);  (b)  Jpo  is  about  6xlCT14  A/cm2  over  a broad  range  of  base 
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width  W (from  15  p to  100  pm)  if  the  back  surface  is  covered  by  ohmic 
contact,  (c)  high  open-circuit  voltage  up  to  700  mV  is  achievable  for 
this  value  of  Jpo.  This  implies  that  we  can  make  thin-base  DF  solar 
cells  to  save  material  and  fabrication  cost  without  sacrificing  the 
conversion  efficiency. 

In  Chapter  3,  we  compared  the  predictions  of  Chapter  2 and  those 
from  other  work  [27]  with  experimental  results  obtained  from  devices 

made  for  this  study.  These  devices  were  non-optimal  in  the  sense  that 
the  DF  region  started  and  ended  within  the  quasi-neutral  base  and 
extended  only  11  pm;  the  thickness  of  the  device  is  38  pm.  Because 
these  devices  showed  reciprocal  slope  m > 2,  indicating  areal 

inhomogeneity  involving  at  the  least  shunt  resistance,  we  estimated  the 
recombination  lifetime  (assumed  independent  of  position  in  the  model) 
from  two  methods  in  which  the  terminals  were  essentially  short 
circuits.  From  the  two  methods  came  the  values:  % = l.o  ps  and 

t = 0.65  ps.  Using  a one-dimensional  model  and  the  non-optimal  doping 
profile  described  above,  we  find  that  despite  these  low  lifetimes  the 
predicted  open  circuit  voltage  is  636  mV  and  625  mV.  These  values 

contrast  with  the  measured  open-circuit  voltage  of  the  best  solar  cell: 

598  mV.  Thus  we  have  two  conclusions.  First,  areal  inhomogeneity  has 

decreased  the  open-circuit  voltage.  Second,  the  results  suggest  the 
immunity  to  the  degrading  influence  of  low  lifetime  of  a DF  cell  design, 
even  for  this  non-optimal  DF  cell. 

Neugroschel  and  Mazer  [53]  found  the  preferential  diffusion  of 
Phosphorus  down  grain  boundaries  in  Wacker  p-polysil icon 
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diffusion  leads  to  n-type  sidewalls  encasing  the  grain-boundary  surface 
and  forming  an  n/p  junction  that  can  extend  into  silicon  a depth  < 18  pm 
[53].  If  this  diffusion  occurs  in  a solar  cell  having  a base  width 
about  20  pm,  the  electrons  inject  into  base  region  from  all  surfaces 
except  the  back  surface  under  forward  bias. 

In  Chapter  4 we  deal  with  the  homogeneous  boundary  condition 
(excess  density  = 0)  at  the  grain  boundary  surfaces.  This  allows  the 
boundary-value  problem  of  Neugroschel  and  Mazer  to  be  treated  as  a 
superposition  of  solutions.  Each  such  solution  is  obtained  from  a 
boundary-value  problem  in  which  the  boundary  density  at  a single  surface 
is  non-zero  whereas  the  boundary  densities  at  all  other  surfaces 
vanish.  This  approach,  slightly  altered  by  using  symmetry,  yields  the 
I-V  characteristics  when  the  recombination  velocity  at  the  back  surface 
becomes  zero,  as  it  would  approximately  if  one  added  an  effective  back- 
surface  field  to  the  structure  proposed  by  Neugroschel  and  Mazer. 

The  use  of  superposition  enabled  by  the  approach  of  Chapter  4 is 
one  of  its  main  advantages.  Other  results  of  Chapter  4 include  the 
demonstration  that  the  computer  simulation  of  Lanza  and  Hovel  [7]  of 
polysilicon  solar  cells  is  questionable  because  the  model  they  use  leads 
to  divergent  current  for  V > 0. 

In  Chapter  5,  we  assess  the  accuracy  of  the  new  method  introduced 
in  Chapter  4.  Following  Gonzalez  [8],  we  solve  the  I-V  characteristics 
of  a grain  having  surface  characterized  by  S = 106  cm/s;  the  result  is 
an  infinite  series,  which  we  truncate  at  600  terms.  The  agreement 
between  a three-term  solution  by  the  method  of  Chapter  4 and  the  600 
term  solution  is  within  a factor  of  2.5  for  the  cases  studied. 
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The  large  discrepancy  between  the  upper-bound  lifetime  in  the  model 
of  Chapter  2 and  the  experimentally  determined  lifetime  of  Chapter  3 
implies  that  defects  other  than  the  fundamental  defects  exist.  We 

recommend  further  work  reducing  the  density  of  these  defects.  We 

recommend  also  an  effort  to  make  devices  having  profiles  more  resembling 
those  of  Chapter  2. 

Both  in  Chapter  4 and  in  Chapter  5,  the  analyses  are  based  on  the 
assumption  of  constant  effective  surface  recombination  velocity  S for 
low-injection  conditions  at  the  edge  of  the  grain-boundary  space  charge 
region  adjacent  to  the  quasi-neutral  region.  In  fact,  S is  position 
dependent  [58].  Theoretical  or  experimental  methods  that  more  precisely 
define  this  position  dependence  are  recommended. 
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